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abstract 


Practical  bounds  on  computer  capacity  and  computation  time 
dictate  a  need  for  simplified  methods  of  airborne  gunnery  ballistics 
computation  particularly  under  conditions  which  involve  large  ballistic 
yaw.  The  projectile  trajectory  data  needed  for  on-board  solution  of 
the  fire  control  problem  can  be  obtained  in  two  ways:  (1)  trajectory 
computation  on  board  the  aircraft,  (2)  curve  fitting  precomputed  tra¬ 
jectory  data.  This  report  contains  a  comparison  of  two  approaches 
to  on-board  trajectory  computation  and  polynomial  curve  fitting  of 
precomputed  trajectory  data. 

A  new,  simplified  set  of  approximate,  large-yaw,  ballistic  tra¬ 
jectory  equations  are  derived  for  airborne  gunnery  applications  in 
which  corrections  for  yaw  drag  and  windage  jump  are  included.  The 
results  of  computations  done  with  the  new  set  of  equations  are  in  ex¬ 
cellent  agreement  with  the  results  of  calculations  done  with  the  approx¬ 
imate  equations  used  at  Eglin  Air  Force  Base  to  generate  trajectory 
tables.  Comparisons  are  made  for  the  20-mm,  M56  round.  A  reduc¬ 
tion  in  computer  time  by  a  factor  of  1  5  to  20  (depending  on  firing  geom¬ 
etry)  over  the  Eglin  set  of  equations  can  be  expected  when  the  simpli¬ 
fied  set  of  equations  is  used.  The  new  set  of  equations  is  almost  as 
simple  as  particle  trajectory  equations  and  is  a  candidate  set  for  on¬ 
board  calculations. 

When  the  yaw-drag  and  windage-jump  corrections  are  used  with 
the  Siacci  method,  a  closed-form  solution  to  the  ballistic  trajectory 
equations  is  obtained.  This  method  is  almost  as  accurate  as  the  simpli¬ 
fied  set  of  equations  above,  except  for  time-of-flight  calculations.  It 
is  sufficiently  accurate  for  air-to-air  applications  out  to  a  range  of 
approximately  3000  ft.  If  time  of  flight  is  not  important,  as,  for  ex¬ 
ample,  in  air-to-ground  applications,  then  this  method  yields  results 
of  sufficient  accuracy  out  to  moderate  ranges  for  most  geometries. 

Also,  this  method  appears  to  be  20  to  1  00  times  as  fast  as  point-mass 
integration  and  500  to  1500  times  as  fast  as  numerical  integration  of 
the  Eglin  approximate  equations  of  motion. 

A  method  of  polynomial  curve  fitting  of  ballistic  lead  angles  for 
air-to-ground  applications  is  presented.  This  method  may  be  useful 
in  situations  where  on-board  calculations  are  not  possible. 
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SECTION  I 


INTRODUCTION 


■s 


1 .  Onboard  Trajectory  Computations  for  Large- Yaw  Situations 

The  use  of  side-firing  or  flexible  gunnery  systems  in  modern  air¬ 
borne  applications  introduces  the  problem  of  large -yaw  ballistics.  Fixed, 
side -firing  guns  are  used  in  the  guns  hip  application  and  the  use  of  flexible 
(turreted)  systems  is  anticipated.  The  problem  also  occurs  to  a  lesser 
extent  in  modern  fixed-gun  fighter  applications  in  high  angle -of-attack 
situations. 

The  problem  of  large -yaw  ballistics  arises  in  airborne  applications 
whenever  the  gun-pointing  direction  lies  significantly  away  from  the  direc¬ 
tion  of  motion  of  the  aircraft.  When  this  occurs,  the  spin  axis  of  the  pro¬ 
jectile,  which  coincides  with  the  gunbore  line,  does  not  coincide  with  the 
direction  of  motion  of  the  projectile.  The  angle  between  the  projectile 
velocity  vector  and  spin  axis  is  known  as  the  projectile  angle  of  attack  or 
the  yaw  angle.  The  effect  of  nonzero  yaw  upon  the  projectile  is  to  pro¬ 
duce  aerodynamic  moments  at  right  angles  to  the  spin  axis.  The  be¬ 
havior  of  the  spinning  projectile  under  the  influence  of  these-  moments  is 
much  like  that  of  a  gyroscope,  or  a  top,  under  the  influence  of  gravity. 

The  projectile  precesses  about  its  velocity  vector. 

Nonzero  yaw  also  induces  aerodynamic  force-S  at  right  angles  to 
the  velocity  vector,  namely  the  lift  and  Magnus  forces.  The  lift  force  is 
in  the  plane  containing  the  projectile  spin  axis  and  velocity  vector  (the 
plane  of  yaw)  and  its  effect  is  analogous  to  ..the  effect  of  the  lift  force  on 
an  airplane.  For  the  case  of  a  spinning  projectile,  however,  the  lift 
force  precesses  with  the  projectile  about  the  velocity  vector.  This  is 
analogous  to  an  airplane  in  a  barrel-roll.  The  Magnus  force  is  perpen¬ 
dicular  to  the  plane  of  yaw  and  behaves  in  a  similar  fashion.  Under  the 
influence  of  these  forces  (mainly  the  lift,  which  is  considerably  larger 
than  the  Magnus  force),  the,  projectile  executes  a  sort  of  spiraling  motion 
as  it  moves  downrange.  For  a  dynamically  stable  projectile,  the  yaw 
angle  rapidly  dies  out  .and  the  spiraling  motion  stops.  The  total  effect  is 
to  deflect  the  tra.jeqtbry  away  from  its  original  direction  by  a  small 
amount.  This  effect  is  the  swerve,  or  the  windage  jump,  and,  as  a  rule 
of  thumb,  it  amounts  to  about  a  milliradian  (mr)  of  deflection  per  degree 
of  projectile- initial  yaw.  For  example,  for  an  F-4  in  a  dogfight,  flying 
at  800  ft/-Sec,  and  with  an  angle  of  attack  of  the  gun  of  15  deg,  the  initial 
yaw  ajrgle  of  the  projectile  is  about  3  deg  and  the  trajectory  will  be  de¬ 
flected  about  3  mr  to  the  right. 

The  aerodynamic  drag  is  also  increased  by  nonzero  yaw.  This 
increase,  referred  to  as  the  yaw  drag,  tends  to  slow  the  projectile  down 
at  a  faster  rate  and  increases  the  time  of  flight.  This  increase  in  time 
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of  flight  is  not  considered  significant  in  air-to-ground  situations,  but  may' 
be  important  in  air-to-air  encounters.  In  air-to-ground  situations,  the 
shape  of  the  trajectory  is  not  altered  to  a  significant  extent  and  it  does  not 
hurt  if  the  projectile  arrives  “a  little  late."  In  air-to-air  situations,  the 
,  target  is  moving  at  a  high  rate  and  small  errors  in  the  time  of  flight  may 
result  in  erroneous  kinematic  predictions  and  large  miss  distances. 

The  motion  of  a  spinning  projectile  with  large  yaw  can,  of  course, 
be  determined  by  means  of  numerical  integration  of  the  six-degree-of- 
freedom  equations  provided  appropriate  aerodynamic  data  is  available. 
However,  the  computation  of  trajectory  tables  may  require  excessive 
computer  time  even  at  ground-based  computer  facilities.  Onboard  solu¬ 
tion  of  the  six-degree-of-freedom  equations  in  real-time  fire-control 
applications  is  not  considered  possible  using  current  technology.  It  may 
be  possible  with  large-scale  integrated  circuits  (LSI),  but  this  problem 
has  not  been  investigated  to  the  author's  knowledge. 

The  trajectory  table  problem  for  the  20-mm,  M56  round  was  solved 
by  use  of  the  set  of  approximate  equations  which  are  documented  in  Ref.  1 
and  which  do  not  require  much  computer  time.  These  approximate  equa¬ 
tions  have  been  written  in  FORTRAN  as  code  R37G  for  the  computer  at 
Eglin  Air  Force  Base  by  Eglin  personnel.  The  equations  are  still  unsuit¬ 
able  for  onboard  use,  however. 

This  report  contains  two  approaches  to  the  problem  of  onboard  tra¬ 
jectory  computation  which  are  reasonably  accurate  for  the  20-mm,  M56 
round  and  which  do  not  require  excessive  computer  time.  These  methods 
involve  approximate  methods  for  treating  windage  jump  and  handling  yaw 
drag. 

In  the  first  method  to  be  described,  a  yaw  drag  correction  is  added 
to  the  point-mass  equations,  and  windage  jump  is  accounted  for  by  utiliza¬ 
tion  of  the  results  of  calculations  done  with  the  equations  of  Ref.  1.  The 
second  method  is  the  Siacci  method  with  a  yaw  drag  correction  and  with 
the  same  windage  jump  correction  used  in  the  first  method. 

The  method  of  treating  windage  jump  was  suggested  for  use  with 
the  Siacci  method  in  Ref.  2,  and  an  investigation  of  its  utilization  is 
reported  herein  and  in  Ref.  3.  A  natural  extension  is  to  use  the  same 
windage  jump  equations  with  point-mass  calculations;  this  was  done  and 
the  results  are  reported  herein  and  in  Ref.  4. 

The  yaw-drag  correction  is  that  of  Sterne  (Refs.  5  and  6).  It  was 
originally  applied  to  the  Siacci  method  and  its  use  has  been  extended  to 
the  modified  point-mass  calculations  herein  and  in  Ref.  4. 

The  equations  of  the  first  method  which  are  derived  herein  and  in 
Ref.  4  are 


2 


where 


.d)[> +  V*kP] 

dt  1 

SP  =  u 


dD  g 

SP  ~  \ 


dQ 
dP  “ 


D 


V  =  u  !  1  -  £D  sin  6  + 


P  is  the  pseudorange,  Q  is  the  gravity  drop,  t  is  the  time,  u  =  dp/dt 
is  the  pseudovelocity,  m  is  the  projectile  mass,  d  is  the  projectile 
diameter,  p  is  the  air  density,  V  is  the  projectile  velocity,  is  the 

velocity  of  sound,  g  is  the  acceleration  due  to  gravity,  Kp  is  the  zero- 

o  — 

yaw  drag,  kQ  and  k  are  constants,  and  9^  is  the  angle  of  elevation  of  P 

above  the  horizontal.  These  equations  are  related  to  a  rectangular  air- 
mass  coordinate  system  by  the  relations 


£  =  P  cos  6  +  C.P 

o  § 


ti  =  P  sin  6  -Q  +  CP 

O  T) 


p 


where  C..P,  C  P,  and  C  P  are  components  of  swerve  (windage  jump); 

»  1  a 

Ct,  C  ,  and  C.  are  calculated  from  initial  conditions  and  are  constants 
S  T1  = 

for  any  given  trajectory. 


To  an  individual  familiar  with  the  six-degree -of-freedom  equations, 
the  simplicity  of  this  set  of  relations  is  obvious.  The  advantage  of  using 
P  as  the  independent  variable  rather  than  t  lies  in  the  fact  that  fewer 
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integration  steps  are  required  to  integrate  out  to  a  given  range,  so  com¬ 
puter  time  is  reduced.  One  might  add  that  these  equations  would  be  part¬ 
icle  trajectory  (point-mass)  equations  if  the  term 

were  deleted.  They  are  applicable  with  greater  accuracy  wherever  part¬ 
icle  equations  are  useful.  The  possibility  of  onboard  numerical  integra¬ 
tion  of  particle  trajectories  has  been  demonstrated  in  the  F-lllD  Mk  II 
avionics  and  in  the  Hot  Line  applications. 

In  the  second  method  (the  modified  Siacci  method)  the  equations  of 
motion  as  described  above  are  replaced  by  the  following  set: 

S{u/aQ)  =  S{uo/aQ)  +  -g-P  + 

a  -  (prf  [m»/s>  -  A4»0/»0i  p] 

>  =  [T<“/s>  - T" Vs*] 

where  C  is  the  ballistic  coefficient,  aQ  is  the  ratio  of  the  speed  of  sound 

at  altitude  to  that  at  sea  level,  <r  is  the  relative  air  density  (<r  =  p/po 

where  pQ  is  the  air  density  at  sea  level),  and  c  is  a  constant.  S,  T,  I, 

and  A  are  functions  tabulated  against  U.  The  term  kQ/2cC  is  the  yaw 

drag  correction  first  derived,  presumably,  by  Sterne.  The  windage  jump 
correction  employed  with  the  Siacci  equations  is  the  same  as  that  de¬ 
scribed  above.  Siacci  ballistic  computations  are  valid  for  large-yaw, 
moderate-range  trajectories  of  the  20-mm,  M56  round.  Applicable  for 
both  air-to-air  and  air-to-ground  fire  control,  the  Siacci  method  yields 
a  substantial  reduction  in  computation  time  over  numerical  integration 
of  the  equations  of  motion.  The  basic  argument  for  implementation  of 
the  Siacci  ballistics  equations  for  onboard  fire  control  computation  is 
that  updating  of  ballistic  lead  angles  may  be  performed  in  a  shorter  time 
interval  than  can  be  accomplished  by  numerical  integration,  thus,  reduc¬ 
ing  the  error  imposed  by  failure  to  obtain  instantaneous  lead  information 
in  a  time -dependent  environment.  In  the  air-to-air  case,  reduction  in 
ballistic  computation  time  affords  more  time  for  processing  of  kinematic 
lead  data  and  other  necessary  functions  of  the  fire  control  computer. 


Siacci  theory  is  an  approximation  to  the  point-mass  equations  of 
motion  of  a  projectile.  The  basic  assumption  is  that  the  velocity  of 
the  projectile  in  the  direction  of  the  initial  velocity  vector  is  equal  to 
the  actual  projectile  velocity  at  all  points  along  the  trajectory.  The 
error  arising  from  this  assumption  is  small  for  short-range,  relatively 
flat  trajectories;  however,  the  error  increases  with  range  and  trajectory 
curvature,  which  limits  the  effective  range  of  Siacci  ballistics.  For  the 
20-ram  round,  fairly  good  computational  accuracy  is  maintained  for  ranges 
approaching  the  effective  limits  of  the  gun. 

2.  Lead  Computation  Modeling 

For  onboard  solution  of  fire  control  problems,  there  are  two  ways 
to  obtain  projectile  trajectory  data.  Trajectories  may  be  numerically 
integrated  on  board  the  aircraft  as  needed  or  else  trajectories  or  bal¬ 
listic  leads  may  be  calculated  beforehand  at  ground  computer  installa¬ 
tions  and  provided  as  an  input  into  the  fire-control  computer  in  the  form 
of  tables  or  curve  fits.  The  method  used  depends  upon  the  type  of  fire- 
control  system,  fire-control  accuracy  requirements,  and  physical  sys¬ 
tem  limitations  such  as  computer  storage  size  and  cycle  time. 

If  approximate  trajectory  calculations  of  one  type  or  another  as 
described  herein  are  considered  to  be  sufficiently  accurate,  onboard  cal¬ 
culations  are  possible  as  demonstrated  in  the  fire-control  system  in  the 
F-lllD  (Mark  II  avaioncs)  and  the  currently  experimental  Hot  Line  sys¬ 
tem.  However,  required  computation  time  is  short.  In  a  system  such 
as  that  of  the  F-lljD,  digital  computation  of  an  entire  trajectory  is 
ideally  instantaneous.  A  practical  bound  on  computation  time,  however, 
may  be  the  reaction  time  of  the  pilot  or  about  one-tenth  of  a  second.  In 
the  Hot  Line  system,  calculations  are  carried  out  in  real  time;  the  avail¬ 
able  calculation  time  is  equated  to  the  time  of  flight.  Multiple  trajectories 
must  be  integrated  at  once,  however.  For  a  3-sec  time  of  flight  and  a 
display  showing  10  "tracer  bullets,"  the  available  calculation  time  is 
0.3  sec  per  trajectory,  if  trajectories  are  computed  serially  on  a  time¬ 
sharing  basis. 

When  onboard  calculations  are  not  feasible,  methods  of  curve  fit¬ 
ting  precomputed  trajectory  data  are  appropriate.  This  investigation  is 
intended  to  aid  in  the  design  of  fire-control  systems  for  flexible  air-to- 
ground  gunnery.  Such  systems  must  have  high  accuracy  to  be  effective, 
particularly  at  relatively  high  operational  altitudes.  Also,  rapid  compu¬ 
tation  of  lead  angles  is  essential  to  this  high  accuracy.  The  computational 
methods  should  also  be  applicable  to  air-to-air  situations. 

The  use  of  curve-fitted  trajectory  data  is  a  tradeoff.  Computer 
time  should  be  minimized,  but  more  core  storage  capacity  maybe  neces¬ 
sary.  The  purpose  of  this  investigation  is  to  determine  some  idea  as 
to  the  amount  of  storage  necessary  and,  in  the  process,  to  discover  the 
best  method  of  handling  the  trajectory  data.  The  best  method  should  use 
the  least  computer  core  space  while  also  using  minimal  computer  time. 
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There  are  a  number  of  schemes  for  handling  ballistics  data  in 
the  airborne  computer  and  a  number  of  methods  for  curve  fitting  it.  The 
most  direct  method  would  be  to  curve  fit  the  trajectories  themselves, 
but  the  quantities  actually  required  in  the  fire-control  problem  are  the 
ballistic  leads  in  one  form  or  another.  One  such  method  is  contained 
herein  (Section  III)  and  in  Ref.  7. 
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SECTION  n 


TRAJECTORY  COMPUTATION 


1.  The  Large-Yaw  Problem 


The  yaw  angle  {angle  of  attack)  of  a  projectile  is  the  angle  between 
its  longitudinal  axis  and  its  direction  of  motion  through  the  air.  When 
a  projectile  leaves  the  muzzle  of  a  gun,  its  longitudinal  axis  orientation 
is  that  of  the  gun  barrel  or  of  the  muzzle  velocity  vector  V  ,  whereas 

the  projectile  velocity  Vo  is  the  vector  sum  of  the  muzzle  velocity  Vm 


and  the  aircraft  velocity  V^;  i.  e. 


V 

o 


As  can  be  seen  from  Fig.  1,  the  angle  between  V  and  V  ,  which  is  the 
yaw  angle  6^,  can  be  large  when  the  angle  between  and  is  large. 

When  the  yaw  angle  is  small,  as  it  usually  is  in  the  fixed  gun  applica¬ 
tion  of  fighter  aircraft  where  and  V  are  nearly  parallel,  the  motion  of 
the  projectile  is  adequately  represented  by  the  simple  equation 


'  eqV  +  mg  (1) 

where  V  is  the  projectile  velocity,  is  the  drag,  m  is  the  projectile 

mass,  and  g  is  the  acceleration  due  to  gravity.  Also 


E  =  p  d  V  Kn  (M) 


(2) 


where  p  is  the  air  density,  d  is  the  projectile  velocity,  Kj^  is  the  zero- 

o 

yaw  drag  function,  and  M  is  the  Mach  number. 


M  =  V/Vsd  (3) 

where  Vgd  is  the  speed  of  sound. 

Equation  {1)  accurately  describes  the  motion  of  a  projectile  mov¬ 
ing  through  the  air  when  the  yaw  is  zero.  Additional  forces  such  as  the 
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yaw  drag,  the  lift  force,  and  the  Magnus  force  result  from  nonzero  yaw, 
so  when  the  yaw  is  too  large,  Eq.  (l)  is  no  longer  a  good  approximation. 


The  lift  and  Magnus  forces  are  at  right  angles  to  the  projectile 
velocity  vector  and  contribute  to  the  swerve  {or  the  windage  jump).  The 
swerve  is  the  deviation  of  the  true  trajectory  from  the  one  described  by 

Eq,  (1).  The  yaw-drag  force  is  in  the  opposite  direction  to  V  and  tends 
to  increase  the  projectile  time  of  flight  between  gun  and  target. 

The  time  history  of  the  yawing  motion  of  a  projectile  is  given  by 
the  equation 


where  H  is  the  angular  momentum  vector  of  the  projectile  and  G  is  the 

aerodynamic  torque.  G  includes  the  Magnus  moment,  the  overturning 
moment,  the  damping  moment,  and  the  spin  deceleration  moment.  It  is 
this  equation  which  presents  the  major  difficulty  in  the  numerical  inte¬ 
gration  of  the  six-degree-of-freedom  equations  because  of  the  gyroscopic, 
nutational, and  precessional  motion  of  the  projectile.  Accurate  descrip¬ 
tion  of  the  nutational  motion  requires  a  very  small  integration  step  size. 

In  order  to  solve  the  large -yaw  problem,  either  the  full  set  of  six- 
degree-of-freedom  equations  must  be  utilized,  or  a  suitable  approxima¬ 
tion  must  be  found.  Such  an  approximation  is  described  in  Ref,  l,  where 
the  nutational  motion  was  eliminated  from  the  equations,  A  further 
simplification  is  documented  here. 

2.  The  Approximate  Equations  Used  in  the  Eglin  Code 


A  yawed,  spinning  projectile  does  not  maintain  a  vertical  plane 
trajectory  as  is  described  by  the  point-mass  equations  of  motion  and  the 
basic  Siacci  equations.  Due  to  the  angular  motion  of  the  projectile,  the 
aerodynamic  forces  cause  it  to  swerve, with  components  in  and  perpen¬ 
dicular  to  the  initial  plane  of  yaw.  Examination  of  calculated  values  of 
the  angular  windage-jump  components  from  the  set  of  approximate  equa¬ 
tions  used  by  Eglin  Air  Force  Base  for  trajectory  tabulation  for  the  20- 
mm,  M56  round  (Ref.  1)  revealed  that  both  angular  components  become 
essentially  constant  after  the  yaw  is  sufficiently  damped.  Further,  both 
components  are  conveniently  approximated  as  being  proportional  to  6q 

with  almost  no  dependence  upon  other  initial  conditions.  This  provided  a 
very  simple  means  of  obtaining  an  accurate  representation  of  fairly  com¬ 
plex  projectile  motion.  The  results  of  Ref.  1  are  summarized  by  the 
following  equations  (see  the  List  of  Symbols). 
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ffiP  =  -  E  P 

(5) 

mQ  =  -  E  Q  +  mg 

(6) 

**  ♦  7  y 

mS  -  -  £5  +  pd  V 
x  x  r 

cos  +  v  Kp  sin  sin  6 

{7> 

mS  =  -  ES  +  pd2V2 

y  y  r 

sin  -  v  Kp  cos  4>'^j  sin  6 

(8) 

cos  6  p  d2V  L  md2  S  ^  .  pd5  „  ?  1  .  2, 

"3t  m  [KL  ~~K  j  KT  +  ^2  [  J  Sln  6 

(9) 

•  _  pd3V2  r  f 

*  ‘  “AW .  [km 

k-  khkt] 

(10) 

where  a  dot  represents  differentiation  with  respect  to  time, 

E  =  pd2VKD{M,6)  (11) 

V  =  J  P2  -ZPQsin  6o  +  Q2  (12) 

and 

4>'  =  <t>  -  4>q  (13) 

Terms  containing  can  be  deleted  from  Eqs.  (9)  and  (10)  since  they 
are  negligible. 

Projectile  coordinates  are  shown  in  Figs.  2  and  3.  In  Fig.  2,  the 
1,2,3  coordinate  system  is  right-handed  with  the  1  axis  along  the  velo¬ 
city  vector  V,  The  1,2  plane  is  a  vertical  plane  containing  V  (V  is  not 
necessarily  horizontal).  The  precession  angie  $  is  a  rotation  about  the 
1-axis  of  2  into  2'  which  is  in  the  plane  of  yaw.  The  plane  of  yaw  is  the 

plane  containing  the  shell  axis  and  the  velocity  vector  V.  The  plane  of 
yaw  rotates  about  V  with  precession  rate  4. 
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Figure  l 

Projectile  Coordinate  System 


jordinates 


rpendicular  to  P 


In  Fig.  3,  f-.rj.f,  is  a  right-handed  coordinate  system.  The  |,r| 
plane  is  vertical  and  contains  the  projectile  initial  velocity  vector  V  . 

The  t,-axis  is  horizontal.  The  pseudorange  P  is  along  V  and  extends 
to  a  point  that  would  be  vertically  above  the  projectile  if  the  swerve  S 
were  zero.  The  gravity  drop  Q  would  be  the  vertical  distance  from  the 
tip  of  P  to  the  projectile  under  the  same  conditions.  Thus,  P  and  Q 
would  describe  the  projectile  trajectory  if  the  swerve  ft  were  zero. 

It  follows  from  Fig.  3  that 

|  =  P  cos  eo  +  (1>K 

■q  =  P  sin  0  -  Q  +  S  (15) 

4  =  (16) 

where  S^,  S  ,  and  arejthe^components  of  S  in  the  system  and 

0  Q  is  the  elevation -angle  of  Vq  or  P  above  the  horizontal.  If  these  rela¬ 
tions  ax-e  Substituted  for  the  components  of  Eq.  (1)  but  with  S  =  0, 

it-iS"seen  that 


«»  • 


mP  =  -E  P 
o 

(17) 

mQ  =  -E^Q  t  mg 

(18) 

The  only  difference  between  these  equations  and  Eqs.  (5)  and  (6)  is  the 
presence  of  E  instead  of  E  in  Eqs.  (5)  and  (6).  But  as  can  be  seen 
from  Eqs.  (2)  and  (11)  ° 

Eq(M)  =  E(M,  6=0)  (19) 

It  follows  that  if  the  time  history  of  6  is  known  from  some  unspecified 
source,  E  could  be  included  instead  of  Eq  and  the  accuracy  of  Eq.  (1) 

could  be  improved.  Also,  if  "3  is  determined  by  some  unspecified  method 
a  complete  solution  to  the  large -yaw  problem  is  available.  It  turns  out 
that  this  is  possible  as  will  be  shown.  It  is  possible  to  effectively  delete 
Eqs.  (7)  through  (10)  from  the  set  and  still  retain  the  accuracy  of  Eqs. 

(5)  through  (13). 


3.  The  Swerve  Approximation- 

The  component^-w  S  are  related  to  S  and  S  by  the  following 
equations  (Ref.^lh'^  X  ^ 


SA  =  -  (S  cos  <J>  -  S  sin  A  )  sin  0 

5  x  ~o  y  ~o  c 

S  =  (S  cos  d>  -  S  sin  A  )  cos  0 

TJ  X  T  o  y  T  o  O 

S.  =  S  sin  A  +  S  cos  A 
^  x  To  y  o 


(20) 

(21) 

(22) 


where  <j>Q  is  the  initial  precession  angle.  Investigation  of  solutions  of 

the  approximate  equations,  i.  e.  ,  Eqs.  (5)  through  (13),  shows  that  the 
windage  jump  components  (in  milliradians ) ,  namely 


103  S^P 

(23) 

103  Sy/p 

(24) 

approach  constant  values  as  the  time  increases,  and  are  essentially 
independent  of  all  initial  conditions  except  initial  yaw  6q.  Furthermore, 

it  is  observed  that  X  and  Y  may  be  approximated  to  sufficient  accuracy 
by  the  relations 


X  =  a6 

o 

Y  =  b  6 

o 

For  the  20 -mm,  M56  round 

a  =  0.  131  mr/deg 
b  =  0.  985  mr/deg 
when  X  and  Y  are  in  milliradians  (mr)  and  6q 


(25) 

(26) 


(27) 

(28) 


is  in  degrees. 


It  follows  that 


where 


£  =  P  cos  0  +  C,P 

o  | 

n  =  P  sin  0  -  Q  +  C  P 

o  T) 


C  =  C^P 


.  3 

-10  6  (a  cos  d>  -  b  sin  A  )  sin  0 

o'  '  o  T  o  o 

10  ^  6  (a  cos  6  -  b  sin  6  )  cos  0 

o'  o 

-  3 

10  6  (a  sin  *  +  b  cos  <J>  ) 

o  ^ o  To 


(29) 

(30) 

(31) 

(32) 

(33) 

(34) 


Thus,  the  swerve  is  given  in  terms  of  trajectory  initial  conditions,  the 
pseudorange,  and  values  of  a  and  b  calculated  from  a  few  selected  tra¬ 
jectories  obtained  by  means  of  the  approximate  equations  of  Ref.  1. 

4.  The  Yaw -Drag  Approximation 

Experience  has  shown  that  the  yaw-dependent  drag  coefficient  may 
be  written  as 


Kd(M,6) 


Kd  (M)  1  +  Kd  2(M)  62 
o  |_  6 


(35) 


There  are  several  ways  to  obtain  the  time  history  of  6,  but  six-degree- 

of-freedom  calculations  show  that  it  is  possible  to  approximate  6  by  a 
relation  of  the  form 


62  =  S2Ke'kP  (36) 

where  K  and  k  are  constants  and  K  =»  1 .  When  yaw  dependence  is  in¬ 
cluded  in  Eq.  (1),  it  is  seen  that 


m 


dV 

-at 


-  E 

o 


1  +  k 

c 


V  +  mg 


(37) 
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where 


k  =  Kn  ,  62  K 
o  DgE  o 


Alternately,  Eqs.  (5)  and  (6)  become 


mQ  = 


**  .lrp  » 

mP  =  -  E  1+ke  P 
o  o 

-E  l  +  k  e“kP  Q 
o  o 


-kP  * 

1+ke  Q  +  mg 


In  the  present  approximation,  these  last  two  equations,  along  with 
Eq.  (12)  for  V,  Eq,  (2)  for  Eq,  and  the  swerve  equations  of  the  last 

section,  replace  the  entire  set  of  approximate  equations  represented  by 
Eqs.  (5)  through  ( 1 3). 

Several  methods  are  available  for  acquiring  values  for  k  and  K. 
Short  of  measuring  5  in  a  spark  range,  the  most  accurate  method  should 
be  to  obtain  yaw  history  data  from  six-degree-of-freedom  calculations. 
Presumably,  Kand  k  do  not  vary  much  from  trajectory  to  trajectory 
(i.  e.  ,  for  different  conditions  of  fire)  so  only  a  few  such  six-degree-of- 
freedom  calculations  need  be  made.  Yaw  data  obtained  in  this  manner 

can  be  plotted  on  semilog  paper  (In  6  4  versus  P);  the  slope  of  the  curve 

is  k  and  the  P  =  0  intercept  is  K6  . 

o 

Short  of  six-degree-of-freedom  calculations ,  Eqs .  (5)  through  (13) 
can  be  used.  Alternately,  the  approximate  equation  for  the  yaw,  Eq.  (9), 
mav  be  written  as 


-  pd2V 


sin  6 


The  damping  term  has  been  deleted  since  it  is  negligible.  Now 


d6  dP 

dp  ~ar 


d6  „  dS 
U-3F  *  VSP 


and  if  sin  6  is  approximated  by  6 ,  it  follows  that 


«  «  .l£-  K 

cTP  m  L 


KT  6 
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Since,  from  Eq.  (36) 


6  =  6  n/k  e-W2 

o 

it  follows  that 

V* .-«■/*  -  -I6 

and  so 

k  “  2  [Kl  *  TST  kt  j  (41) 

K  can  be  taken  equal  to  one,  or  the  value  to  be  defined  below  by  Eq.  (42). 

Another  approach  is  to  borrow  the  results  of  the  Siacci  theory  as 
formulated  by  Sterne,  (Refs.  5  and  6)  for  which  the  average  squared  yaw 

T 

6  is  given  by 


,  *°  '  'll  62  e-2«P 

s  -  1  o 


(42) 


where  represents  the  static  stability  factor  which  is  given  by 


s  = 
o 


a2n2 

4Bpdi^KM 


(43) 


and 


c'  + 


s  -  1 

o 


(44) 


where 


W  +  T-  kh] 


(45) 


17 


and 


tstV 


B  is  the  transverse  moment  of  inertia  of  the  projectile  and 


Pq  =  0.  076474  lb/ft“ 


is  the  air  density  at  sea  level. 

5.  The  Modified  Point-Mass  Method 


With 


and  by  means  of 


u 


=  P 


d 

dt 


(.  \  •  • 

_Q\=  PQ  -  PQ 

P  /  P ~ 

it  can  be  shown  that  the  complete  set  of  equations  to  be  solved  is 
mi  =  -Pd2VKD^V/Vad)  Tl  +  koe-kPja 


P  =  u 


and 


±/_Q.\  =  _S_ 

dt  \  u  J  u 


•here 


V  = 


V 


PZ  -  2PQ  sin  e  +  QZ 
o 


(46) 


(47) 

(48) 


(49) 


(50) 


18 


In  many  applications  it  has  been  found  convenient  to  make  the  independent 
variable  P  rather  than  t.  This  is  done  as  follows 


du  du  dP  du 

sr  ~  3P"  ar* =  uhf 

and  hence,  Eq.  (47)  becomes 

m  "HP  *  -l*d2vKDo(V/V»d)[‘+''<)e'kP 


From  Eq.  (48) 


dt  1 
SF  ~  u 


and  from  Eq.  (49) 


But 


Q  _  Q  _  dQ/dt  .  dQ 
u  p  dP/dt  “  HF" 

and  so 


If  we  define 


D 


dQ 

cfF- 


it  follows  that 


dD 

■3F 


(51) 


(52) 


(53) 


(54) 


(55) 
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Also 


Q 


dQ  dP 

=  -dP  "HT  " 


u 


dQ 

SP 


and  Eq,  (50)  becomes 

V  =  -  2D  sin  0q  +  D2  (56) 

The  new  set  of  equations  to  be  solved  numerically  are  Eqs.  (51),  (52), 
(54),  (55),  and  (56). 

6.  The  Siacci  Equations  and  Correction  Terms 

b.  1  Basic  Siacci  Equations  -  The  basic  Siacci  equations  are  also 
approximations  to  point-mass  equations  of  motion.  The  method  provides 
a  closed-form  solution  to  projectile  ballistics  after  the  tabulation  of 
numerically  integrated  functions  of  the  zero-yaw  drag  coefficient.  The 
derivation  of  the  equations  has  been  well  documented  since  their  con¬ 
ception  in  the  1880's  by  F.  Siacci  of  Italy  (Refs.  2,  6,  and  8).  No  effort 
will  be  expended  here  in  rederivation. 

The  tabulated  Siacci  functions  necessary  for  solution  of  the  Siacci 
trajectory  equations  are  obtained  by  numerical  integration  of  the  follow¬ 
ing  differential  equations. 


dS 

"an  = 

•Gim 

(57) 

dT 

dtf  = 

-  uG(CT 

(58) 

dl 

-an* = 

-  V— 

U  G(U) 

(59) 

dA 

I(U) 

(60) 

"air  _ 

-  wr 

U  is  given  by  U  =  u/aQ  where  u  is  the  velocity  along  the  initial  velocity 
vector  and  aQ  is  the  ratio  of  the  speed  of  sound  at  the  firing  altitude  to 
that  at  sea  level.  =1116.45  ft/sec  is  the  speed  of  sound  at  sea  level 
for  the  standard  atmosphere.  is  the  zero-yaw  drag  coefficient  and 

2  ° 

g  =  32. 174  ft/sec  is  the  acceleration  due  to  gravity.  The  functions  S, 

T,  I,  and  A  are  tabulated  for  incremental  values  of  U.  The  differential 
equations  are  integrated  from  an  arbitrary  upper  bound  of  U  to  an  arbi¬ 
trary  lower  bound.  A  tabulation  of  these  values  from  U  =  7,  000  ft/sec 
to  U  =  500  ft/sec  for  the  20-mm,  M56  round  may  be  found  in  Ref.  9. 

The  Siacci  tables  are  independent  of  altitude,  as  has  been  shown  in 
Refs.  2,  6,  and  8. 

The  trajectory  equations  are  given  as  follows 

P  =  £-  [S  (u/ao)  .  S(uo/ao)j  ((>2) 

Q  '  (^)  [A(“/a°)  *  A(“°/ao)  *  I(“o/a„)-&  P]  <63> 

*  *  -^r  [T(u/ao)  -T(v/ao)]  <64> 

where  t  is  the  time  of  flight  of  the  projectile.  P  is  directed  along  the 

initial  velocity  vector  uq  =  VQ,  Q  points  vertically  down,  and  u  is  defined 

as  u  =  dP/dt.  The  basic  Siacci  approximation  to  the  point-mass  equations 
of  motion  is  that 


V  =  u 

As  is  evident,  this  approximation  does  not  hold  for  extremely  large 
values  of  Q.  The  symbol  C  is  the  ballistic  coefficient  and  it  is  the  rel¬ 
ative  air  density  at  the  firing  altitude. 


C 


m  2 

- - — k -  l  lb/m. 

144  a  X 


(65) 


where  m  is  the  mass  of  the  round  in  pounds,  and  d  is  the  projectile 
diameter  in  feet, 


6.  2  The  Yaw-Drag  Correction  -  Since  the  Siacci  equations  were 
derived  from  point- mats s  equations,  the  projectile  trajectory  is  assumed 
to  be  in  a  vertical  plane  with  no  effect  of  yaw  on  the  projectile  drag. 
These  assumptions  are  reasonably  accurate  for  very  small  initial  yaw 
angles,  but  corrections  must  be  introduced  for  large-yaw  ballistics. 

According  to  the  theory  of  Sterne,  the  yaw  drag  is  taken  into 

"T 

account  by  the  approximation  that  the  average  squared  yaw,  6  ,  is  an 
inverse  exponential  function  of  P  given  by 


7?  o  ^  -Z<r  cP 

o  =0  — - - —  e 

o  s  -  1 
o 


(66) 


as  has  been  previously  explained.  The  correction  term  for  P,  according 
to  Sterne,  is  given  by 


AP  -  -  Kd  2  \  62  dP 
6  0 


jr 

I 


(67) 


where  6  is  the  yaw  angle  and  Kn  ,  is  the  yaw-drag  coefficient.  Combin¬ 

es2 

ing  Eqs.  (66)  and  (67)  with  6^  substituted  for  yields 


AP  =  - 


2<rc 


[  -2acP  "I 

>-e  J 


(68) 


where 


s _  -  1/2  , 

k  =  — - r~  Kn  ,  6 

o  s  -  1  Df2  o 

o  0 


(69) 


The  equation  for  S(u/aQ)  may  be  obtained  by  adding  the  correction  term, 

Eq.  (68),  to  the  right-hand  side  of  Eq.  (62)  and  rearranging  terms.  In 

most  cases,  the  yaw  damps  rapidly  enough  for  e  ^<rc^>  to  be  neglected 
for  practical  ranges.  Therefore 

k 

S(u/a  )  =  S(u  [ a  )  +  -1-P  +  (70) 

o  0  0  C  2c  C 


22 


Thus,  it  is  apparent  that  if  the  initial  conditions  u  and  6  are  known 

o  o  * 

along  with  the  atmospheric  quantities,  ballistic  parameters,  and  P,  then 
Q  and  t  are  solved  by  examining  the  tabulated  values  of  T,  1,  and  A. 

6.  3  Curve  Fit  of  the  Siacci  Functions  -  Storage  of  an  extensive 
listing  of  Siacci  tables  is  cumbersome  and  requires  a  large  core-storage 
capacity.  For  this  reason,  the  Siacci  functions  of  the  20-mm,  M56  round 
were  curve  fitted  as  fourth-degree  polynomials.  The  tables  were  broken 
up  into  sections,  and  each  section  was  curve  fitted  at  equally  spaced 
values  of  the  independent  variable*  Reasonably  good  curve  fit  accuracy 
was  obtained.  The  tables  that  were  curve  fitted  are  contained  in  Ref.  10. 
The  curve  fitting  was  done  prior  to  the  generation  of  the  new  tables  in 
Ref.  9,  but  the  two  sets  of  Siacci  tables  are  essentially  identical. 

The  tabulated  values  of  S  were  curve  fitted  as  a  function  of  U, 
while  the  values  of  T,  I,  and  A  were  fitted  as  functions  of  S.  The  in¬ 
dependent  variable  U  extended  over  the  range  of  7,000  to  1,  000  ft/sec. 
Though  some  of  the  trajectories  computed  for  the  comparison  in  Section 
V  have  a  terminal  velocity  slightly  below  1,000  ft/sec,  extrapolation  of 
the  curve  fits  does  not  appear  to  induce  severe  error  for  the  ranges 
considered.  Caution  should  be  exercised  for  ranges  in  excess  of  5,  000  ft. 
For  longer  ranges  and  improved  accuracy,  a  more  accurate  curve  fit  of 
the  tabulated  values  of  Ref.  9  would  be  appropriate.  The  coefficients  of 
the  curve  fits  used  in  this  analysis  are  given  in  Appendix  I. 

7.  Initial  Conditions 


A  comparison  has  been  made  between  calculations  done  with  the 
simplified  methods  described  herein  and  those  done  with  the  approximate 
equations  of  Ref.  1.  Because  of  the  interest  in  air-to-ground  fire-control 
problems  at  the  time  this  work  was  done,  air-to-ground  geometry  was 
adopted  for  this  comparison.  The  equations  for  the  necessary  initial- 
condition  computations  and  coordinate  transfoi  mations  are  given  in  this 
section.  The  results  of  these  calculations  are  given  in  Section  IV. 

It  is  assumed  that  an  aircraft  flies  straight  and  level  at  a  given 
altitude  over  a  sea-level  target  range  and  that  the  aircraft  carries  a  gun 
turret  armed  with  an  M61,  20-mm  cannon.  The  ballistics  for  this  round 
are  well  known  {Ref.  11).  The  turret  gimbaling  is  assumed  to  be  of  the 
depression,  traverse  type  shown  in  Fig.  4.  The  x,  y  plane  is  horizontal 
and  z  is  vertical.  The  aircraft  flies  along  the  x-axis,  and  y  points  out 
the  left  wing.  The  depression  angle  E  is  a  positive  rotation  about  the  y- 
axis  and  rotates  x  into  x'.  The  muzzle  velocity  V  (i.e.  ,  the  gun)  is  in 

the  x',y  plane,  and  A  is  the  angle  between  x'  and  V  .  Initial  conditions 

for  computations  are  defined  by  the  equations  below  along  with  the  trans¬ 
formation  from  the  £ ,  r\,  t,  system  to  the  x,  y,  z  system. 
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Components  of  VQ  =  along  x,y,  and  z  are 


V 

V. 

+  V  cos  A 

cos  E 

ox 

A 

m 

V 

V 

sin  A 

oy 

m 

V 

-  V 

cos  A  sin 

E 

oz 

m 

Also 


(71) 

(72) 

(73) 


V  =-t/vf  +  2V.  V  cosAcosE  +  V2 
o  yj  A  A  m  m 


(74) 


and 


6  =  arccos 

o 


V  +  VA  cos  A  cos  E 
m  A 

- v - 


(75) 


The  elevation  angle  0  of  V  (shown  positive  in  Fig.  5)  is  defined  by  the 
equations  °  J 


sin  0  = 

o  V 

o 


(76) 


and 


cos  0 


o  V 


JiiL 


(77) 


where 


The  angle  B  (in  Fig.  5)  is  defined  by 


V 


xy 


(78) 


(79) 
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cos  B 


(80) 


I 


ox 

xy 


The  coordinate  system  is  defined  as  shown  in  Fig.  5.  The  £ 

axis  lies  along  the  projection  of  V  in  the  horizontal  x,y  plane,  r)  is  coin¬ 
cident  with  z,  and  £  completes  the  right-handed  set. 


The  initial  value  of  the  precession  angle  4>  can  be  found  by  examina¬ 
tion  of  Figs.  2  and  5.  At  time  t  =  0,  the  1,2,3  coordinate  axes  of  Fig,  2 

lie,  respectively,  on  the  Vq,  z',  and  r,  axes  of  Fig.  5.  Then  4>q  is  the 

angle  between  z'  and  the  projection  of  V  in  the  z',{L  plane.  It  follows 
that  m 


cos  <t> 
To 

sin  cf> 


(cos  A  cos  E  cos  B  +  sin  A  sin  B)  sine 
_ _ _ o 

sin  6 

o 

cos  A  cos  E  sin  B  -  sin  A  cos  B 
sin  8 

o 


+  cos  A  sin£  cos 


(81) 

(82) 


The  coordinates  of  a  point  in  transform  into 


x  =  %  cos  B  +  £  sin  B  (83) 

y  =  £  sin  B  -  £  cos  B  (84) 

z  =  *1  •  (85) 

in  the  x,  y,  z  system. 

In  the  calculations,  Eq,  (42)  was  used  to  define  the  projectile  yaw 
history.  Values  used  as  ballistic  parameters  for  the  20mm,  M56  round 
are 


c'  =  0.0023  ft"1 
c"  =  0.000105  ft’1 


13.  2 
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m  =  0.224156  1b 
d  a  0.  06562  ft 

and 


where 

s  =  3.5655 
s 

Also 


V  =  3300  ft/sec 
m  ' 

Values  of  A  and  B  used  are 

A  =  0.  00013  lb-ft2 
B  =  0.00096  lb-ft2 

Also 


N  =  12,  1349. 6  rad/sec 

Worst-case  calculations  were  also  done.  The  worst  case  occurs  at 
high  airspeeds  when  the  gun  is  positioned  so  the  initial  yaw  is  maximum. 
As  seen  from  Fig.  6,  maximum  yaw  is  given  by 


E'  =  6  +  90° 

o 

It  is  noted  that  the  modified  point-mass  calculations  and  the  modified 
Siacci  calculations  were  done  at  different  times  and  the  coordinate  systems 
used  in  the  two  calculations  were  different.  Only  the  coordinate  system 
used  for  the  modified  point-mass  calculations  is  reported  herein  6ince  the 
two  systems  are  equivalent  and  do  not  affect  the  results. 
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SECTION  m 


POLYNOMIAL  CURVE  FITTING  OF  BALLISTIC  LEADS 


The  method  presented  herein  for  curve  fitting  trajectory  tables  has 
been  used  in  the  past  for  bombing  trajectories  (Ref.  12).  Extension  of 
the  method  to  the  curve  fitting  of  firing  tables  is  an  academic  matter, 
but  must  be  tested  to  show  that  a  high  degree  of  accuracy  may  be  ob¬ 
tained  over  a  fairly  wide  range  of  coverage  for  an  airborne  turreted  gun 
system.  Bombing  tables  may  generally  be  represented  by  one  dependent 
variable,  such  as  the  trail  or  the  spherical  earth  range-to-target,  tabu¬ 
lated  as  a  function  of  altitude,  true  airspeed,  etc.  Due  to  the  necessity 
for  a  much  greater  delivery  accuracy  for  gunfire,  it  will  be  necessary  to 
curve  fit  at  least  two  dependent  variables,  i.e. ,  ballistic  lead  angles 
relative  to  the  line  of  sight  (LOS)  to  the  target.  For  kinematic  lead  pre¬ 
diction  and/or  wind  corrections,  it  may  also  be  necessary  to  curve  fit 
the  time  of  flight  of  the  projectile.  The  latter  variable,  however,  will 
not  be  needed  for  stationary  ground  targets  in  still  air,  and  it  may  be  pos¬ 
sible  to  approximate  the  time  of  flight  to  sufficient  accuracy  for  kine¬ 
matic  lead  prediction  and/or  wind  corrections  by  means  other  than  a 
polynomial  curve  fit. 

The  matrix  method  presented  in  the  following  sections  will  be 
studied  from  the  standpoint  of  the  air-to-ground  situation.  The  target 
will  be  assumed  to  be  stationary  at  sea  level,  the  aircraft  will  be  assum¬ 
ed  to  be  flying  straight  and  level,  and  it  is  assumed  crosswind  is  neg¬ 
ligible. 

Four  independent  variables  are  used  to  curve  fit  the  depression 
and  deflection  ballistic  lead  angles  relative  to  the  LOS  to  the  target. 

They  are  relative  air  density,  aircraft  velocity,  depression  angle  of  the 
LOS,  and  deflection  angle  of  the  LOS.  Aircraft  altitude  is  extracted 
from  the  relative  air  density.  Thus,  the  range  to  the  target  may  be  de¬ 
rived  from  the  altitude  and  LOS  angles.  Two  curve  fits  are  necessary 
to  obtain  both  the  depression  and  deflection  ballistic  lead  angles. 

Equally  spaced  values  of  the  independent  variables  are  used  for 
programming  convenience.  However,  the  derivation  of  the  matrix  equa¬ 
tions  does  not  contain  this  restriction. 

The  Siacci  method  (Refs,  2,  6,  and  8)  was  used  to  generate  the 
data  to  be  curve  fitted.  Although  this  method  is  not  as  accurate  as  the 
numerical  integration  of  the  projectile  equations  of  motion,  it  is  suffi¬ 
ciently  accurate  for  current  exploratory  computations  and  is  much  more 
economical.  The  same  type  of  curve  fits  can  be  obtained  for  the  more 
accurate  ballistic  data  when  they  are  needed. 
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1.  Multidimensional  Matrix  Derivation  for  Exact  Fit  to  the  Data  Points 


Consider  one  dependent  variable  b  and  four  independent  variables 
w,  x,  y,  and  z.  The  values  of  b  are  known  fqr  I,  J,  K,  and  L  different 
values  of  w,  x,  y,  and  z,  respectively.  Then,  b  may  be  expressed  as  a 
polynomial  function  of  the  independent  variables  and  constant  coefficients 

Aijki  aS  shown  in  E<P 


b 


I-1,J-14K-1.L-1 

h 

i,  j,k,f 


i  j  k  l 
w  xJy  z 


(86) 


The  values  of  must  be  found  by  use  of  the  values  of  b  correspond¬ 

ing  to  the  various  combinations  of  w,  x,  y,  and  z.  Therefore 


pqrs 


1-1,  J-l.K-l.L-l 

l 

i.  j,  k,  f 


\  a  i  j  k  i 

Li  ljkl  p  q y  r  s 


(87) 


where 


P  =  0,  1,  2,  1-1 

q  =  0,  1,  2, - ,  J-l 

r  =  0,  l,  2,  K-l 

s  =  0,  1,  2, - ,  L-l 


In  matrix  form,  W  -  (W.  }  .  etc.  .  whe 

\  ip/' 


re 


W.  =  w  ,  etc. 

ip  p' 


Define  the  inverse  of  the  independent  variable  matrices  as  follows 

C  =  W"1,  D  =  X*1,  E  =  Y'1,  F  =  Z"1 
so  that 
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V 

L 

p 


W.  C  =6. 
ip  pa  xa 


y*.  Db 

L  jq  qb 
q 

y  y  e 

l_j  kr  re 


(88) 


\  Zfs  Fsd  "  6|d 


Thus 


^  b  C  D  E  F  , 
pqrs  pa  qb  re  sd 

P.  q.  r,  s 


=  )  A. .  6.  6  ,6,  6,  , 

ljki  la  jb  kc  id 

i,  j.k.i 


(89) 


The  6  terms  are  elements  of  the  identity  matrices;  hence 


A  ,  .  =  Y  b  W'1  X'i  Y-i  Z~\  (90) 

abed  Pqrs  pa  qb  ec  sd  ' 

P.  q.  r.  s 

Equation  (90)  is  the  symbolic  form  of  the  general  coefficient  matrix  and 
may  be  easily  expanded  to  N-dimensions  by  inspection.  Substitution  of 
Eq,  (90)  into  Eq.  (86)  along  with  specific  values  of  the  independent  variables 
yields  the  desired  dependent  variable. 

Normalizing  the  independent  variable  matrices  will  simplify  the 
matrix  algebra  and  reduce  round-off  error  acquired  in  matrix  inversion. 
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Suppose  the  independent  variables  are  <r  (relative  air  density),  (air¬ 
craft  velocity),  A ^  (LOS  deflection  to  target)  and  E^  (LOS  depression  to 
target).  Act,  AV^,  AAj,,  and  AE,p  are  equal  intervals  between  succes¬ 
sive  variations  in  the  independent  variables.  The  subscript  m  indicates 
the  value  of  the  variable  at  the  mid- point  of  the  region  covered  by  each 
variable.  The  normalized  variables  used  in  the  curve  fit  are  then  de¬ 
fined  as  follows 


w 


cr  -  tr 


m 


A<r 


x  = 


V  -  V 
A  v  Am 

ATT 


y  = 


z  = 


AT  *  ATm 
ET  -  ETm 

~ae: 


Tm 


The  matrices  W,  X,  Y,  and  Z  will  be  of  the  form 


(91) 


W 


1  .  .  . 

•  1 

.  .  . 

*  wI-l 

2 

2 

.  .  . 

*  WI-1 

• 

• 

• 

• 

• 

1-1 

1-1  .  . 

1 

*  WI-1 

(92) 


The  order  of  matrix  (92)  is  equal  to  the  number  of  data  points  and  will  de¬ 
pend  upon  the  degree  of  the  polynomial.  The  polynomial  b(w,x,y,  z)  is 
of  degree  1-1  in  w.  Hence,  I  data  points  are  required  and  the  order  of 
W  is  I.  The  W  matrix  is  shown  in  Eqs.  (93)  through  (96)  for  1=2,  3,  4, 
and  5. 
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For  I  =  2 


For  1=3 


Fori  =  4 


W  = 


and  for  1=5 


W  = 


1 

-0.  5 


1 

0.  5 


W  = 


1  1 
0  1 
0  1 


- 

1 

1 

1 

1 

-1.5 

-0.5 

0.  5 

1.5 

2.25 

0.25 

0.  25 

2.  25 

-3.  375 

-0.125 

0.  125 

3.  375 

w  = 


1 

-2 

4 

-8 

16 


1 

-1 

1 

-1 


1 

0 

0 

0 

0 


1 

1 

1 

1 

1 


1 

2 

4 

8 

16 


(93) 


(94) 


(95) 


(96) 
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It  may  also  be  necessary  to  normalize  the  dependent  variable  to  reduce 
errors  associated  with  round  off.  If  is  the  value  of  the  dependent 

variable  a  when  the  independent  variables  are  <r  ,  V_  ,  A_  ,  and 

m  um  I  m 

E„  ,  and  N  is  a  normalizing  constant  which  may  be  chosen  to  reduce 
1  m  a 

the  values  of  the  dependent  variables  to  near  unity,  the  dependent 
variable  to  be  curve  fitted  is  taken  to  be 

a  -  a 

b  =  -  (97) 

a 

If  the  values  of  a  are  not  excessively  large  or  small,  N  =  1  will  be 
ade  quate . 

2.  Programming  and  Storage  for  the  N-Dimensional  Arrays 

Since  both  b  and  A  are  multidimensional  arrays,  it  is  convenient 
to  store  their  values  in  a  linear  array  fashion  (our  CDC  3200  computer 
cannot  handle  four -dimensional  arrays  with  the  currently  used  FORTRAN 
IV  system).  The  following  sample  program  should  serve  to  relax  pro¬ 
gramming  complexity  which  occurs  when  the  multidimensional  equations 
are  broken  up  into  two-dimensional  matrix  algebra. 

It  is  assumed  that  the  variables  have  been  normalized,  and  the 
values  of  the  dependent  variable  at  equal  intervals  are  given.  The  coef¬ 
ficient  matrix  can  be  obtained  by  use  of  the  following  program. 

PROGRAM  CURVE  FIT 

C  4 -DIMENSIONAL  CURVE  FIT 

C  INDEPENDENT  VARIABLES  XI,  X2.X3.X4 

C  DEPENDENT  VARIABLE  Z 

C  COEFFICIENT  MATRIX  A 

C  NUMBER  OF  DATA  POINTS  FOR  EACH  VARIABLE 
C  RESPECTIVELY,  I,  J,  K,  L 

DIMENSION  Z(625),  Xl(5,5),  X2(5,  5).  X3(5,5),  X4(5,5),  A(625) 
READ  10,  I,  J,  K,  L 
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10  FORMAT  (411) 


READ  20,  {(X1(N,M),  M  =  1,1),  N  =  1,1) 

READ  20,  ((X2(N,  M),  M  =  i,J),  N  =  1,J) 

READ  20,  ({X3(N,  M),  M  =  1,K),  N  =  1,K) 

READ  20,  {(X4(N,M),  M  =  1,1,),  N  *  i,L) 

20  FORMAT  (1  OF  8.1) 

KL  =  1*J*L*K 

READ  20,  (Z(N),  N  =  1,  KL) 

C  SUBROUTINE  INVERT  (X)  INVERTS  X.  THE  X  ARRAY  STORAGE 

C  IS  USED  FOR  STORAGE  OF  THE  INVERSE 
CALL  INVERT  (XI) 

CALL  INVERT  (X2) 

CALL  INVERT  (X3) 

CALL  INVERT  (X4) 

DO  40  II  =  1,1 
DO  40  JJ  =  1,  J 
DO  40  KK  =  1,  K 
DO  40  LL  =  1 ,  L 

N  =  (II-l)* J*K*L  +  (JJ-1)*K  *L+  (KK-1)*L  +  LL 
A(N)  =  0. 

DO  30  IF  =  1,1 
DO  30  IQ  =  1,  J 
DO  30  IR  =  1,K 
DO  30  IS  =  1,L 
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M  =  (IP-1)  *J  *K*L,  +  (IQ-1)  *K*  L.  +  (IR-1) *  L  +  IS 
A(N)  =  A(N}  +  Z(M)  *  X1(IP,II)  *  X2  (IQ,  JJ)  *X3(IR,KK) 
*  X4  (IS, LiL) 

30  CONTINUE 
40  CONTINUE 

PUNCH  50,  (A(N),  Na  1,  KL) 

50  FORMAT  (4E20.  12) 

END 


It  should  be  noted  that  the  X  matrices  should  be  read  in  by  varying 
columns  and  then  rows.  The  values  of  the  dependent  variable  Z 

ljkf 

should  be  read  in  while  varying  the  parameters  in  reverse  order,  i.  e. , 

I ,  k,  j,  and  i.  The  output  of  matrix  A  will,  therefore,  be  in  the  same  * 
order  as  Z  was  input.  Although  the  program  was  set  up  for  4-dimen¬ 
sional  curve  fit  with  a  maximum  of  five  data  points  for  each  independent 
variable,  it  may  be  easily  expanded  to  a  higher-degree  fit  by  insertion 
of  the  required  parameters  and  an  increase  in  the  size  of  the  dimension¬ 
ed  arrays. 

3.  Evaluation  of  the  Accuracy  of  the  Curve  Fit  for  Air-to-Ground  Fire 
of  a  20-mm  Gun  "  “ 


To  determine  the  range  through  which  the  independent  variables 
may  be  varied  and  to  what  degree  they  must  be  curve  fit,  it  is  necessary 
to  estimate  the  degree  of  each  variable  and  test  the  curve  fits  at  points 
between  the  mesh  of  data  points.  This  leads  to  a  trial-and-error  se¬ 
quence  in  which  an  effort  is  made  to  optimize  the  curve  fit  for  a  pre- 
specified  accuracy.  Any  desired  accuracy  may  be  obtained  by  decreas¬ 
ing  the  range  of  coverage  and/or  increasing  the  degree  of  the  curve  fit. 
The  problem  then,  is  one  of  attempting  to  minimize  the  number  of  curve- 
fit  coefficients  necessary  to  describe  all  of  the  ballistic  lead  angles 
within  the  effective  coverage  of  the  gun  system. 

A  turreted,  20-mm  gun  was  selected  for  a  quantitative  evaluation  of 
the  curve -fit  method  because  of  the  availability  of  the  ballistic  data 
necessary  for  trajectory  calculations.  The  Siacci  method  was  used  to 
determine  the  lead  angles  to  be  curve  fit,  rather  than  numerical  integra¬ 
tion  of  the  equations  of  motion,  to  reduce  the  computation  time  required. 
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Coordinate  systems  are  defined  as  follows:  In  Figs.  7  and  8  the 
x,  y,  z  coordinate  system  is  earth  fixed  with  x  and  z  horizontal  and  y 
vertical.  The  angle  B  is  rotated  clockwise  about  the  y-axis,  and  is 
measured  between  the  x-axis  and  the  |-axis.  The  £-axis  lies  along  the 

projection  of  the  vector  uo  in  the  horizontal  plane.  The  r^-axis  is  co¬ 
incident  with  the  y-axis  and  the  £-axis  completes  a  right-hand  £  Sy8. 
tem.  JT he  angle  9q  is  measured  between  uq  and  the  £-axis.  The  vectors 
Vm»  Va>  and  UQ  are»  respectively,  the  projectile  muzzle  velocity,  the 
aircraft  velocity,  and  the  projectile  initial  velocity,  ‘  where 


In  Fig.  7,  Eq  is  the  depression  angle  of  the  gun  measured  between 
the  x-axis  and  the  projection  of  Vm  on  the  x,  y  plane.  A^,  the  deflection 
angle  of  the  gun,  is  the  angle  between  and  the  x,y  plane. 

In  Fig.  8,  Az  and  E^  are  the  azimuth  and  elevation  angles  of  the 

gun,  respectively.  A^  is  measured  between  the  projection  of  V  in  the 

x,z  plane  and  the  x-axis,  and  E_  is  measured  between  V  and  the  x  z 

u  m  * 

plane . 


gun 


First,  the  azimuth -elevation  system  (Fig.  8)  was  used  for  the 
pointing  angles,  but  some  difficulty  arose  in  the  curve  fitting  of  the 
azimuth  lead  angle.  A  fourth-order  curve  fit  in  all  four  variables  pro¬ 
duced  errors  that  were  generally  less  than  1  mr,  but  were  as  high  as  10mr 
at  some  points  on  the  curve  fit.  The  relative  air  density  was  varied  from 
0.  98  to  0.  86  (approximately  700  to  5000  ft  in  altitude)  while  the  aircraft 
velocity  azimuth  LOS  angle,  and  elevation  LOS  angle,  respectively  were 
varied;  (200  to  500  ft/sec),  (0  to  120  deg,  and  (-45  to  -75  deg). 


In  an  effort  to  decrease  the  maximum  error  of  the  azimuth  curve 
“fj  gun  lead  angles  were  changed  to  the  depression-deflection  system 
(Fig.  7).  This  angular  arrangement  proved  to  be  considerably  easier  to 
fit  and  corresponds  to  a  more  natural  gimbal  system  for  a  turreted  air- 
to-ground  gun.  The  angular  coverage  in  depression  angle  Eq,  and  de¬ 
flection  angle  AG,  are,  respectively,  -30  to  -135  deg  and  -30  to  30  deg. 

The  aircraft  velocity  ranged  between  200  and  800  ft/sec,  and  the  relative 
air  density  between  0.98  and  0.  86. 


Figures  9,  10,  11,  and  12  represent  plots  of  the  depression  and  de¬ 
flection  lead  angles  vs  target  LOS  angles  as  calculated  by  the  Siacci 
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Figure  7 

Depression  -  Deflection  System 
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Figure  10 

Depression  Lead  Angle 
vs  LOS  Deflection  to  Target 
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Figure  11 

Deflection  Lead  Angle  vs  LOS 
Depression  to  Target 
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Figure  12 

Deflection  Lead  Angle  vs  LOS 
Deflection  to  Target 


method.  These  plots  were  used  to  aid  in  the  first  estimate  of  the  degree 
of  the  polynomial  required  for  the  two  independent  variables  A^,  (LOS 

deflection  to  target)  and  E^.  (LOS  depression  to  target).  Figure  9  sug¬ 
gests  that  AE  (depression  lead)  may  be  approximated  by  a  fifth-degree 
polynomial  in  E^,.  The  accuracy  of  the  curve  fit  is  arbitrary,  but 

directly  affects  the  number  of  coefficients  required,  A  convenient  tol¬ 
erance  of  ±2  mr  was  selected,  well  within  the  gun  pointing  accuracy  of 
current  airborne  fire-control  systems.  The  original  computer  code  was 
set  up  for  computation  of  the  coefficients  for  a  polynomial  of  no  higher 
than  fourth-degree  in  any  of  the  independent  variables.  After  several 
unsuccessful  attempts  to  obtain  fourth-degree  curve-fit  accuracy  with  a 
±2  mr  tolerance  for  AE  over  the  entire  range  of  E^,,  A^,,  cr,  and  , 

the  curve  fit  was  separated  into  two  regions.  The  variable  E^,  ranged 

between  30  and  90  deg  in  the  first  region,  and  between  90  and  135  deg  in 
the  second  region.  The  remaining  three  variables  extended  over  their 
entire  range  in  both  regions  of  the  curve  fit.  After  examination  of  Figs. 

10,  11,  and  12,  and  some  trial- and-error  experimentation,  it  was  found 
that  a  good  curve  fit  could  be  obtained  with  variables  of  the  following 

degrees:  AE  {  a(2),  VA(2),  AT(3),  ET(4)}  over  two  regions,  and 

AA{<r(2),  vaU).  E^,(3)}.  The  variable  AA  required  only  one  curve 

fit.  The  number  of  data  points  required  for  each  independent  variable  is 
equal  to  the  degree  plus  one.  Therefore,  each  section  of  AE  required 
3  X  3  X  4  X  5  =  180  points,  and  AA  required  3X3X4X4  =  144  points 
for  a  total  of  2  X  180  +  144  =  504  points.  This  ia  also  the  total/number  of 
coefficients  in  the  polynomial  curve  fits. 

4 .  Comparison  of  the  Siacci  Curve  Fit  With  20-mm,  Large -Yaw  Tra¬ 
jectory  Tables 

This  section  contains  a  comparison  of  the  Siacci  method  of  ballistic 
computation  with  the  approximate  equations  of  motion  of  Ref,  1,  (Also, 
see  Section  IV  of  this  report.)  The  approximate  equations  were  used  in 
the  computation  of  M56,  20-mm  trajectory  tables  compiled  by  the  Air 
Force  Armament  Laboratory  (Ref.  13).  The  purpose  of  the  comparison 
is  not  to  show  that  the  Siacci  method  is  a  highly  accurate  means  of  per¬ 
forming  ballistic  computations  for  large -yaw,  long-range  trajectories, 
but  rather,  to  show  that  the  curve  fit  of  the  Siacci  data  could  be  used  for 
limited  angular  coverage  and  limited  range  for  onboard  ballistic  lead 
computation.  A  more  realistic  approach,  however,  is  to  curve  fit  the 
lead  angles  obtained  from  the  approximate  equations  of  motion  or  any 
other  method  which  is  better  suited  to  large -yaw  ballistic  computations. 
This  is  not  done  here  because  of  the  relatively  lengthy  time  required  for 
ballistic  lead  computation  by  methods  requiring  small  step-size,  numer¬ 
ical  integration. 
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Selected  trajectories  in  the  firing  tables  (Ref.  13}  were  used  for 
comparison  with  the  Siacci  curve  fit.  The  three  components  of  range  to 
the  target  in  the  tables  were  used  to  resolve  the  LOS  angles.  The  leads 
were  then  computed  by  subtracting  the  gun  angles  from  the  LOS  angles. 
Contrary  to  the  firing  tables,  the  depression  angle  for  the  curve  fit  is 
considered  positive  when  below  the  horizontal. 

The  table  in  Appendix  HI  contains  a  complete  list  of  the  depression 
and  deflection  lead  angles  as  obtained  from  both  the  Siacci  curve  fit  and 
the  firing  tables.  Comparison  was  made  between  the  lead  angles  for 
gun  depression  angles  of  45,  65,  85,  and  105  deg.  Reasonable  agree¬ 
ment  was  obtained  for  the  smaller  depression  angles,  but  for  105  deg 
there  was  a  large  deviation  between  Siacci  calculations  and  the  firing 
tables.  This  may  be  due  to  round-off  error  in  the  firing  tables  at  near 
90  deg  depression  angles  or  computational  errors  in  the  direction- 
cosines  for  an  initial  velocity  vector  depression  angle  greater  than 
90  deg.  Only  data  for  the  smaller  depression  angles  will  be  considered. 
Figures  13  and  14  are  plots  of  the  difference  between  Siacci  calculations 
and  the  firing  tables.  It  should  be  noted  that  a  biased  deviation  between 
the  two  methods  may  be  essentially  eliminated  by  adjustment  of  the 
ballistic  coefficient  used  in  the  Siacci  calculations,  or  by  adjusting  the 
first  coefficient  of  the  curve  fit.  By  simple  inspection  of  the  graphs,  it 
appears  that  an  average  deviation  for  both  the  depression  and  deflection 
lead  angles  is  approximately  ±0.  5  deg  (±8.  7  mr)  with  the  bias  removed. 

The  results  of  this  section  should  be  compared  with  those  of 
Section  IV  where  Siacci  computations  are  also  compared  with  results  of 
calculations  using  the  approximate  equations  of  motion  of  Ref,  1.  In 
both  instances,  the  same  set  of  approximate  equations  was  used.  In  this 
section,  the  difference  between  leads  calculated  by  use  of  the  two  methods 
may  be  as  much  as  1  deg, whereas  the  difference  indicated  in  Section  IV 
is  usually  less  than  1  mr.  The  suspected  reason  for  this  discrepancy  will 
be  discussed  in  Section  IV. 

5.  Variables  not  Included  in  the  4-D  Curve  Fit 


The  four -dimensional  curve  fit  presented  here  is  complete  only  if 
we  assume  straight-and-level  flight,  a  sea-level  target,  and  wind  and 
atmospheric  deviations  from  standard  are  negligible.  The  flight  altitude 
and  target  altitude  may  be  accurately  taken  into  account  by  increasing 
the  number  of  independent  variables  in  the  curve  fit  or  by  curve  fitting 
for  several  target  altitudes  and  interpolating  for  points  between  curve 
fits.  Alternately,  consider  the  following  argument.  Near  sea  level,  the 
air  density  p  is  approximately  given  by 


P(z) 


P  e 
ro 


-hy 
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ERROR(AA)  -  deg 


Figure  14 

Difference  Between  the  Deflection  Lead  Angle 
Curve  Fit  and  the  Firing  Tables 
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where  is  sea-level  air  density,  y  is  altitude  above  sea  level,  and  h  is 
a  constant.  If  y^.  is  target  altitude  and  y^  is  aircraft  altitude,  the  rela¬ 
tive  air  density  between  the  aircraft  and  the  target  is 


*hyA 

P(yA>  =  Ppe  =  e‘h(yA  '  yT> 

pTy^T  “  -hyT 

p  e 
*o 


But,  this  is  the  same  as 


p*yA  “  yT*  -MyA-yx) 


Use  of 

P(yA) 

^  "  P^T7” 

in  the  curve  fits  should  be  accurate  enough  for  most  purposes  since  p(y) 
does  not  change  much  over  a  few  hundred  feet.  The  curve  fits  given 
here  are  for  targets  at  or  near  sea  level.  Curve  fits  for  combat  zones 
at  higher  altitudes  can  easily  be  obtained. 

For  most  air-to-ground  applications,  the  only  parameter  needed 
in  addition  to  the  ones  listed  above  would  be  the  dive  angle.  As  neces¬ 
sary,  any  other  attitude  parameters  could  be  treated  by  means  of  ap¬ 
propriate  coordinate  transformations.  Effects  of  wind  and  atmospheric 
deviations  can  probably  be  handled  by  some  form  of  differential  cor¬ 
rections  to  lead  angles.  In  most  cases,  however,  these  effects  will  be 
small  and,  possibly,  negligible. 


SECTION  IV 


SUMMARY  AND  CONCLUSIONS 


1.  Modified  Point -.Mass  Method 


A  new,  simplified  set  of  approximate,  large-yaw,  ballistic  tra¬ 
jectory  equations  has  been  developed  which  are  in  excellent  agreement 
with  the  set  of  approximate  equations  used  at  Eglin  Air  Force  Base  to 
generate  trajectory  tables.  When  the  new,  simplified  set  of  equations  is 
used,  a  reduction  in  computer  time  by  a  factor  (depending  on  firing 
geometry) of  15  to  20  over  the  old  set  can  be  expected.  This  new  set  of 
equations  is  almost  as  simple  as  the  particle  trajectory  equations  used 
in  the  F-lllB  Mk  II  avionics  system  package  and  is  a  candidate  set  for 
airborne  calculations  in  some  applications. 

A  comparison  has  been  made  between  computations  done  with  the 
simplified  method  described  herein  and  those  done  with  the  approximate 
equations  used  at  Eglin  AFB  (Ref.  1).  The  results  of  these  calculations 
are  summarized  in  Tables  I,  II,  and  III  for  three  aircraft  altitudes, 

2000  ft,  5000  ft,  and  8000  ft,  respectively.  Calculations  were  done  for 
three  values  of  aircraft  velocity  (200  knots,  300  knots,  and  400  knots) 

and  two  values  of  gun  traverse  angle  A  (0°  and  30°).  Three  values  of  gun 

depression  angle  E  (30°,  60°,  and  90°)  appear  in  Tables  I  and  II  and  two 

values  (60°  and  90°)  appear  in  Table  III.  The  initial  yaw  angle  6q,  the 

projectile  time  of  flight,  and  the  range  from  the  gun  to  the  target  are 
tabulated  for  each  trajectory. 

Coordinates  £  and  £  represent  the  projectile  impact  point  on  the 
ground  (at  sea  level)  as  calculated  with  the  old  approximate  method,  and 
and  A£  represent  the  amount  by  which  calculations  done  with  the  new 
simplified  method  differ  from  those  done  with  the  old  method  (old  minus  sim¬ 
plified).  As  seenby  the  pilot  (Fig.  1 5),£  is  measured  downrange  on  the  ground 

along  the  projection  of  Vq  (the  initial  projectile  velocity  vector)  in  the 

horizontal  plane,  and  £  is  measured  to  the  right.  The  radial  difference 
is  the  square  root  of  the  sum  of  the  squares  of  and  Agreement 

between  the  two  methods  is  excellent.  It  is  doubtful  that  the  old  approxi¬ 
mate  calculational  method  compares  with  actual  measured  trajectories  as 
well  as  these  two  calculations  agree. 

Worst-case  calculations  were  also  done.  The  worst  case  occurs 
for  high  aircraft  velocities  when  the  gun  is  positioned  so  the  initial  pro¬ 
jectile  yaw  (angle  of  attack)  is  maximum.  Maximum  yaw  occurs  when 

Vq  =  V m  +  V is  perpendicular  to  V^,  e.g.  ,  for  straight-down  fire. 

Computations  were  done  for  altitudes  of  2000  ft,  5000  ft,  and  8000  ft  with 
an  airspeed  of  600  knots  (just  under  Mach  1),  a  traverse  angle  A  of  Odeg, 


50 


Figure  1 5 

Air-to-Ground  Coordinates 


Comparison  Between  the  Simplified  Approximate  Theory  and  the  Old  Approximate 
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and  a  depression  angle  of  107,9  deg.  In  these  situations  the  initial  yaw  is 
17,  9  deg.  The  results  of  these  computations  are  given  in  Table  IV.  The 
biggest  component  of  the  radial  difference  is  seen  to  be  The  amount 

of  this  difference  for  all  three  altitudes  is  about  3  mr,  and  this  can  be 
corrected  in  the  simplified  theory  by  adjusting  the  swerve  parameters. 

In  all  instances,  the  a;rcraft  is  assumed  to  be  flying  straight  and 
level.  The  firing  geometry  is  that  of  Fig.  15. 

The  complete  new  set  of  differential  equations  for  trajectory  com¬ 
putation  is  listed  here  for  convenience  as  follows: 


where 


m  =  -Pd2vKDo(V/Vsd)[1+ko''1'P] 

dt  l 

HF  =  u 

dD  g 

"3F“  "  “T 

u 


dQ 

"3F 


D 


V 


u  k! 
\ 


1  -  ZD  sin  0  +  D 


Definitions  of  symbols  are  in  the  List  of  Symbols  at  the  front  of  this 
report.  It  is  noted  that  these  equations  would  be  particle  trajectory 
equations  if  the  term 


1 


k  e 
o 


-kP 


were  deleted. 

Integration  of  these  equations  yields  gravity  drop  Q,  time  of  flight 
t.  and  projectile  velocity  V  versus  pseudorange  P.  Projectile  position 
in  air  mass  coordinates  is  gi/en  by 

R  =  Ft  Q  +  S 
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where  S  is  the  swerve.  It  is  verified  herein  that  S  can  be  approximated 
adequately  by  the  relation 


S  =  CP 


where  C  is  given  in  terms  of  initial  conditions  and  two  parameters  a  and 
b.  Components  of  Rin  the  air  mass  system  are 


£  =  P  cos  0  +  C..P 

o  4 

ti  =  P  sin  0  +  C  P  -  Q 

1  O  T) 


4 


where 

Ce  =  -10'3  6  /a  cos  d>  -  b  sin  d>  \  sin  G 
4  o\  ^o  /  o 


C  =  10  3  6  /a  cos  A  -  b  sin  d>  \  cos  G 
T)  o  \  To  /  o 


■a 

C,  =  10  6  /a  sin  A  +  b  cos  <{>  \ 

4  o \  To  o / 

■»  a 

The  4-axis  lies  along  the  horizontal  projection  of  P,  q  is  vertical  positive 
up,  and  4  completes  a  right-handed  coordinate  system.  Angles  Bq,  6q, 

and  4>q  are  respectively  the  elevation  angle  of  P  above  the  horizontal, 

the  initial  yaw  angle,  and  the  initial  precession  angle. 

The  swerve  approximation  is  based  upon  the  observation  that  wind¬ 
age  jump  parameters  X  and  Y  (Refs.  1 ,  2,  and  1 1)  calculated  by  use  of  the 
Eglin  code  R370  are  approximately  constant  after  the  projectile  has  traveled 
about  1000  ft, and  are  related  to  the  yaw  by  approximately  X  =  a6  and 

Y  =  b5  .  ° 

o 

Values  of  k  and  kQ  in  the  differential  equation  for  u  can  be  obtained 

in  several  different  ways  as  is  discussed  in  the  main  body  of  this  report. 
The  method  used  for  the  calculations  recorded  in  this  section  is  borrowed 
from  the  old  theory  of  Sterne  (Refs,  5  and  6).  The  equations  are  as  follows: 
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o 


■J 

t 

M 


a2n2 


k  = 


c 1  = 


s  -  1/2 
o 


Kt 


s  -  1  D-2  o 

o  o 


»  d2  r 

M* 


Tm“  I  +  KH 


] 


P  d 

'o  K 

~~2m~ 


c  =  c'  + 


a  -  1 

o 


k  =  2<rc 


Trajectory  initial  conditions  are  calculated  from  the  following 
relations  which  are  derived  in  Section  II. 


V  = 

V  .  +  V  cos  A  cos  E 

ox 

A  m 

V 

V  sin  A 

oy 

m 

V  = 

-  V  cos  A  sin  E 

oz 

m 

V  =  /'v.2  +  2V.  V  cos  A  cos  E  +  V  2 

o  A  Am  m 

V  +  V.  cos  A  cos  E 
t  „  m  A 

0  =  arccos  — - — — — 

°  o 

V  =  /v"2  +  V z 

xy  V  ox  ov 


|BiS£?m:£-£  ■ 


s 


sin  6  = 


V 

oz 

V 


cos  6  = 

o 

-x.y 

V 

o 

sin  B  = 

V 

— 21. 
v 

xy 

cos  B  = 

V 

ox 

V 

cos  <j> 


xy 

{cos  A  cos  E  cos  B  +  sin  A  sin  B)  sin  0  +  cos  A  sin  E  cos  6 

_ _  _ o _ o_ 

sin  6 


cos  A  cos  E  sm  B  -  sin  A  cos  B 

sm  <b  =  - < — . . . — — - 

sm  6 


In  addition  u  =  V  ,  t  =  0,  D  =  0,  and  Q  =  0  when  P  =  0, 
o 

The  relation  between  f,rj,4  coordinates  andx,y,z  coordinates  is 


x  =  §  cos  B  +  £,  sin  B 
y  =  £  sin  B  -  £  cos  B 
z  =  ri 

The  x,  y ,  z  system  is  right  handed,  with  x  along  the  flight  path  and  z  vertical. 
2.  Modified  Siacci  Method 


The  modified  Siacci  method  is  applicable  to  both  air-to-air  and  air- 
to-ground  airborne  fire-control  computation.  The  basic  difference  be¬ 
tween  the  two  modes  of  operation  is  the  kinematic  lead  prediction  require¬ 
ment  for  a  moving  target,  occurring  chiefly  in  the  air-to-air  encounter. 
Accurate  kinematic  prediction  requires  an  accurate  estimate  of  projectile 
time  of  flight  as  vrell  as  an  almost  instantaneous  solution  to  the  fire-con¬ 
trol  problem.  The  Siacci  solutions  provide  a  decided  advantage  in  com¬ 
putation  time  over  numerical  integration  of  even  the  simplest  forms  of  the 
equations  of  motion.  It  does  not,  however,  yield  equivalent  accuracy  in 
the  time-of-flight  calculation.  For  the  20-mm,  M56  round,  the  Siacci 
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error  in  time  of  flight  is  less  than  0.01  sec  for  trajectory  ranges  up 
to  3,000  ft,  but  increases  to  values  greater  than  0.  05  sec  for  ranges 
beyond  5,  000  ft.  The  time  of  flight  is  influenced  by  variations  in  air 
density,  thus  increasing  the  error  resulting  from  the  Siacci  method  for 
relatively  high  altitude,  air-to-ground  fire  due  to  failure  to  account  for 
variations  from  the  firing  altitude  air  density.  This  effect  is  not  of 
consequence  when  altitude  variations  are  not  excessive  (air-to-air),  and 
the  time  of  flight  is  not  usually  needed  to  a  high  degree  of  accuracy  for 
air-to-ground  fire  control.  Corrections  for  the  air-density  variation 
have  been  derived  (Ref.  2)  but  will  not  be  considered  here. 

For  the  purpose  of  comparison,  air-to-ground  trajectories  for  a 
turreted  gun  were  computed  by  the  modified  Siacci  method  and  the  ap¬ 
proximate  equations  of  motion.  Air-to-ground  trajectories  were  selected 
over  air-to-air  because  of  the  larger  projectile  yaw  angles  which  are  a 
major  shortcoming  of  the  basic  Siacci  theory.  The  intention  is  to  dem¬ 
onstrate  that  the  modified  Siacci  equations  are  capable  of  accurate  large- 
yaw  ballistic  calculations  for  moderate  ranges.  Calculations  were  made 
for  a  set  of  initial  conditions  similar  to  those  used  for  the  modified  point- 
mass  equations  computations,  and  the  results  are  to  be  found  in  Tables  V, 
VI,  and  VII.  As  before,  altitudes  of  2000  ft,  5000  ft,  and  8000  ft,  and 
aircraft  speeds  of  200  knots,  300  knots,  and  400  knots  are  used.  Gun 
depression  angles  are  30  deg,  60  deg,  and  90  deg  and  traverse  angles 
are  0  deg  and  30  deg.  The  gun-pointing  error  is  seen  to  be  less  than 
1  mr  for  all  ranges  below  5000  ft,  no  greater  than  2  mr  for  ranges  up 
to  8000  ft,  and  as  high  as  5.  1  mr  for  ranges  slightly  above  9000  ft.  The 
largest  errors  occur  at  the  higher  altitudes  and  longer  ranges  due,  in 
part,  to  the  variation  in  air  density  and  the  breakdown  of  the  approxima¬ 
tion  V  =  u  (Section  II).  However,  the  modified  Siacci  method  yields  a 
highly  accurate  representation  of  the  approximate -equation  solution  for 
yaw  angles  as  high  as  11.  56  deg  and  some  ranges  greater  than  8000  ft. 

The  results  of  worst-case  calculations  are  shown  in  Table  VIII.  The 
results  are  about  the  same  as  those  for  the  worst- case,  modified  point- 
mass  calculations  of  Table  IV,  and  the  reason  is  the  same;  the  windage 
jump  approximation  can  be  improved  for  these  conditions. 

As  a  last,  but  important,  word  on  comparison  of  the  modified 
Siacci  method  with  more  accurate  forms  of  ballistic  equations,  three 
FORTRAN  ballistics  computer  codes  were  compared  for  relative  speed 
of  computation.  The  modified  Siacci  code  was  compared  with  modified 
point-mass  integration  and  the  approximate  equations  of  motion.  Though 
no  attempt  was  made  to  optimize  the  coding  of  the  three  programs,  it 
is  believed  that  ratios  of  the  computation  times  should  be  fairly  repre¬ 
sentative  of  what  would  be  experienced  in  an  airborne  fire-control  com¬ 
putation  .  The  relative  times  vary  from  one  trajectory  to  another  as 
functions  of  range  and  step  size  of  the  numerical  integration;  but,  for 
the  trajectories  considered,  the  Siacci  method  appears  to  be  20  to  100 
times  as  fast  as  the  modified  point-mass  integration,  and  500  to  1500 
times  as  fast  as  numerical  integration  of  the  approximate  equations  of 
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motion.  Representative  calculation  times  obtained  by  use  of  the  CDC 
3200  computer  are  0.  01  sec  for  the  Siacci  method,  0.  5  sec  for  point- 
mass  integration,  and  10  sec  for  the  approximate  equations  of  motion. 

If  the  modified  Siacci  approach  is  to  be  used  in  airborne  calcula¬ 
tions,  the  same  equation  set  as  that  listed  above  for  the  modified  point- 
mass  equations  may  be  used,  except  the  set  of  differential  equations  is 
deleted  and  replaced  with  the  Siacci  equations  as  follows: 


k 

S(u/a  )  =  S(u  /a  )  +  -~-  P  +  — 

'  o  o'  o  C  2  c  C 


•  (*)‘  [ 
|t{u 


A(u/aQ)  -  A(uQ/ao) 


- 

-  T.  P 


C 
<r  a. 


/aQ)  -  T(uo/ao) 


] 


The  functions  S,  T,  A,  and  1  are  precomputed  at  a  ground-based  facility 
and  are  made  available  onboard  in  tabular  form  (Refs.  2  and  9)  or  as 
curve  fits  (Appendix  I). 


3.  Polynomial  Curve  Fitting  of  Ballistic  Lead  Angles 

Polynomial  curve  fitting  of  ballistic  leads  by  exact  fit  to  the  data 
points  appears  to  have  promise  in  some  instances  for  turreted  gun 
systems  that  have  large  fields -of-coverage  and  require  fast  computation 
speed.  The  reduction  of  the  amount  of  data  that  must  be  stored  if  curve 
fitting  is  used  in  place  of  tables,  and  the  computation  time  reduction  if 
it  is  used  in  place  of  numerical  integration,  are  obvious  advantages. 

The  curve  fits  of  ballistic  leads  obtained  in  this  study  by  means  of  the 
Siacci  method  indicate  the  ease  of  accurate  curve  fitting.  The  curve - 
fit  coefficients  are  listed  in  Tables  XVIII,  XIX,  and  XX  of  Appendix  IV. 
Extension  of  the  method  to  the  curve  fitting  of  ballistic  leads  obtained 
by  more  exact  and  complex  forms  of  trajectory  computation  should  pre¬ 
sent  no  problem. 

The  accuracy  of  the  curve  fits  obtained  in  this  study  was  tested  by 
computation  of  the  lead  angles  at  all  points  midway  between  all  data 
points  of  the  grid  used  to  generate  the  coefficient.  Comparison  of  the 
leads  derived  from  Siacci  calculations  with  the  curve  fit  gives  an  in¬ 
dication  of  the  curve -fit  accuracy  obtained.  The  complete  results  of 
this  comparison  are  listed  in  Tables  XIV,  XV,  and  XXI  of  Appendix  II, 

A  summary  of  the  error  analysis  is  given  in  Table  IX  below. 
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Table  IX 


Accuracy  of  Curve  Fits  as  Compared  with  Siacci  Calculations 


Lead 

Angle 

Depression 

Angle  Interval 

Maximum 
Error  (mr) 

Average 
Error  (mr) 

A  E 

30°  s  ET  £  90° 

1.  50 

0.  61 

90°  £  ET  <  135° 

1.28 

0.41 

AA 

30°£ET£l35° 

0.  89 

0.  35 

0.  86  £  <r  £  0.  98  200  ft/sec  £  £  800  ft/sec 

130°£  At  £  30° 


The  accuracy  could  be  increased  or  relaxed  to  any  desired  level 
which  would  show  a  corresponding  increase  or  reduction  of  the  number 
of  coefficients  required  for  the  curve  fit.  With  the  accuracy  obtained 
in  this  example,  360  coefficients  are  required  for  AE  and  144  coeffi¬ 
cients  are  required  for  AA.  Thus,  a  total  of  504  coefficients  would 
have  tc  be  stored  in  the  onboard  computer  of  the  fire-control  system. 

In  the  case  of  tabular  trajectory-table  storage,  many  times  this  amount 
of  data  would  need  to  be  stored  to  obtain  equivalent  accuracy. 


Sensitivity  of  Ballistics  Calculations 


In  Section  III,  a  comparison  was  made  between  the  results  of  cal¬ 
culations  performed  using  the  Siacci  method  and  computations  made  by 
use  of  the  approximate  equations  of  Ref.  1.  A  similar  comparison  was 
made  in  this  section,  and  the  results  of  Section  III  were  poor  relative  to 
these.  The  reason  for  this  discrepancy  is  thought  to  be  due  to  the  fact 
that  this  laboratory  did  not  use  the  same  numerical  values  for  certain  bal¬ 
listic  parameters  (such  as  the  projectile  mass)  as  did  Eglin,  since  the 
approximate  calculations  of  Section  III  were  performed  at  Eglin  Air  Force 
Base,  whereas  those  of  this  section  were  done  at  this  laboratory.  This 
indicates  that  the  results  of  ballistics  calculations  are  sensitive  to  small 
errors  in  ballistics  parameters,  and  one  might  expect  differences  of 
the  order  of  one  degree  between  calculations  made  by  different  labora¬ 
tories.  Note  should  be  taken  of  the  implications  of  this  statement  with 
regard  to  fire-control  system  design. 


APPENDIX  I 


CURVE -FIT  COEFFICIENTS  FOR  THE  SIACCI  FUNCTIONS 


The  Siacci  functions  (S,  T,  I,  and  A)  were  curve  fitted  as  fourth- 
degree  polynomials  over  segments  of  the  data  listed  in  Ref.  10.  The 
tabulated  values  of  S  were  curve  fitted  as  a  function  of  U,  while  the 
values  of  T,  I,  and  A  were  fitted  as  functions  of  S  as  indicated  be¬ 
low. 


S  =  a  +  a,X  +  a,X2  +  a„X3  +  aJC4 
o  1  Z  3  4 


[  T,  I,  or  A]  =  b  +  b,Y  +  b,Y2  +  b,Y3  +  b„Y4 
*  O  1  L  3  4 


The  functional  forms  of  the  independent  variables  X  and  Y  are  given 
by 

/  U  -  U  \ 

x  =  4(irm?x-Tr,r  ) 

\  max  mm  / 

and 

*  *  •) 


The  terms  U  ,  U  .  .  S  ,  and  S  represent  the  upper  and  lower 
max  min  max  mm 

bounds  of  the  curve  fit  intervals  as  shown  in  the  tables.  The  values  of 
a^  and  b^  are  also  given  in  Tables  X  through  XIII  for  specified  intervals 

of  U  and  S. 
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Coefficients  for 
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APPENDIX  II 


ACCURACY  OF  THE  BALLISTIC  LEAD  POLYNOMIAL 

CURVE  FITS 


The  following  tables  are  the  results  of  a  comparison  of  the  poly¬ 
nomial  curve  fits  of  the  depression  and  deflection  ballistic  lead  angles 
with  the  values  derived  by  Siacci  calculations.  The  comparison  is  made 
at  the  mid-point  between  all  data  points  of  the  grid  used  to  generate  the 
coefficients.  The  average  error  should,  in  reality,  be  less  than  the 
average  error  shown  here  because  the  curve  fit  is,  in  general,  more 
accurate  near  the  data  points.  The  error  columns  in  the  tables  are  the 
Siacci  values  minus  the  curve -fit  values  in  milliradians. 
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Deflection  Lead  Angle  Comparison 
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APPENDIX  ID 


COMPARISON  OF  THE  SLACCI  CALCULATIONS  WITH  THE 
EGLIN  20-ram  FIRING  TABLES 

Selected  values  of  the  independent  variables  were  used  to  com¬ 
pare  the  curve  fit  of  the  Siacci  calculations  and  the  Eglin  firing  tables 
for  the  M56,  20-ram  round.  From  the  table,  an  evaluation  of  the  ac¬ 
curacy  of  the  Siacci  calculations  was  made. 
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Table  XVII 

Comparison  of  Lead  Angles  Derived  from  the  Firing  Tables  and  Siacci  Calculations 
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APPENDIX  IV 


BALLISTIC  LEAD  POLYNOMIAL  CURVE -FIT  COEFFICIENTS 


The  following  tables  contain  the  coefficients  obtained  from  the 
direct  polynomial  curve  fit  of  the  ballistic  lead  angles  derived  from 
Siacci  calculations  for  a  turreted,  air-to-ground,  M61  gun.  The  target 
is  assumed  to  be  at  sea  level  and  the  aircraft  is  flying  straight  and 
level.  Effects  due  to  aircraft  attitude  variation  and  wind  velocity  are 
not  included.  The  four  independent  variables  cover  the  following  ranges: 

200  ft/sec  ^  VA  <  800  ft/sec  -30  deg  ^  A^,  2S  30  deg 

0.86  £  ir  :£  0.  98  30  deg  E^  S  135  deg 


Tables  XVIII  and  XIX  list  the  depression  leads  for  line-of-sight  de 
pression  angles  from  30  deg  to  90  deg  and  90  deg  to  135  deg,  res¬ 
pectively.  Table  XX  contains  the  deflection  leads  for  all  values  of  the 
independent  variables  within  the  ranges  specified  above.  The  coef¬ 
ficients  A^.^  are  defined  by  the  following  polynomial. 


Lead  Angle  =  Anii+AiH2 


0  0  0  1 
w  x  y  z 


+  Ajj22  w^x^y^z*  +  .  .  . 


+  A. 


KLMN 


K-l  L-l  M-l  N-l 
w  x  y  z 


where  w,  x,  y,  and  z  are  defined  by  Eq.  (91)  in  Section  III.  The  values 
of  the  subscripts  (K,  L,  M,  N)  are  equal  to  the  number  of  data  points  for 
each  independent  variable  used  in  the  curve  fit. 

The  coefficients  are  presented  in  standard  FORTRAN  "E"  format 
As  an  example,  5.  01326E  -  04  is  equivalent  to  5.  01326  X  10 


Table  XVIII 

Depression  Dead  Coefficients  for  Depression  Angles 
Between  30  deg  and  90  deg 


1  -7.2970E  00  00  2.8360E-01  -i.7?81E-f>2  9.0373E-03 

2  4.2043E-0?  3-004?lr>02  1.1567E-02  -8.8609E-03  -S.1010E-03 

3  -4.  i*oSi--02  4.i921fc-<)2  -1.49<M)p-G2  4.73Q9F-03 

4  -6.5454E-03  -l.blllE-02  -7.1654E-03  3.5201E-03  2.26S9E-G3 

1  -?.9S0?sE-01  1  •  ‘*87Sh— t»l  -7.4276E-02  *>.  1 103F-G2  -1  ,84«6E-n2 

2  7. 764SE-0S  -4*309(|f -OS  8.62ftSfc-t)6  -3.2416E-05  l.*8?HE»05 

3  -8*13b?F-ft2  (<«()HcSSt--02  ~3.23u3fc'-0?  I.3I94F-02  -?.S4n?F-03 

4  -1.4067E-04  7.S98SF-0S  I.0bM2E-0,3  4,14S?E-06  -3.1043E-06 

1  l.l407fc”0|  -y.SSS7F-.42  4,9Sf*7h-0?  -3.1909E-0?  «*4?7?F-03 

2  -3.2371E-0S  2.0329E-05  -7,Bti31fc-06  1.74J3E-05  -7.9922E-06 

3  4. 4 1  l <-t  —  *>>'  -S.0  74OL-02  2.0n3nfc-i)2  8.Hb30F-04  -2.37 4?F-r)3 

4  #>.2  701 E-0S  -4.0344E-0S  -2.4880E-07  -7.3372E-07  S.4«99E-07 


(a)  j  =  1 
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Table  XVIII  (continued) 


i 


1 

2 

3 

4 

5 

4.*»327E  00  -f>.9487E-0) 

1.4?f,7E-0l 

1.  inmfc-o?  -2.44go*;-n3 

5.1910E-02 

3.6395E-02 

1.3816E-02 

-U1003E-02  -f5.2*5h‘ir-0* 

4.1fl71F-01  -4*5511E-«1 

2.38S9E-02 

-2.*>l22t-;)  \ 

1  .'54  1  ^  -04 

9.4147E-03  -1.9396E-0? 

-R.6224E~03 

4.477St-0  3 

2.40  73E-0  1 

1.17SQE-03  -1.8987E-0? 

1.1704E-02 

8.8681E-05  -4*4527E-0«i 

-3.8147E-0S 

-3.307^p-0S 

2.  1  h.hit  -0*5 

6.7529P-03 

2.2411E-03 

-1.9402E-03 

1.|?V'*IM  -J  .  1  ‘)4*r -0  i 

1.SH17E-04 

1.2337E-04 

1.0526E-05 

3*  39i)3t  ”0*1  -4.97«mt-ob 

1.3749F-03 

1.25M2E-0? 

-4.5S43E-03 

-4.  i 

;  -jw  -o  .1 

3.6737E-0S 

1.9512E-0S 

-3.9166E-06 

2.0124E-OS  -1.0417^-0*5 

2.7037E-03  -1.8087E-03 

2.4717F-03 

-2.  17141  - >  < 

*5 . 4  i  4  /l  —  i)4 

7.06S6E-05  -6.3654E-05 

3.2791E-08 

-6. 7H97K-0 ' 

4.4  79<5E-<)7 

<b)  j  =  2 
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Table  XVIH  (concluded) 


! 

1 

2 

3 

4  5 

1  -1  .S9/0E-02 

-1.3622E-02 

3.1076E-04 

3.4909E-Q4  R.1994E-0S 

2  1 .6427E-02 

1.1141E-02 

4.1125E-03 

-3.5204F-03  -1.9573E-03 

3  -2.H492E-02 

I.S45/E-02 

4.12S1E-03 

-9.0144E-04  — 1 . 064QE  —04 

4  -3.7644E-03 

-6.5608E-03 

-P.5935E-03 

1.4933E-03  8.9395E-04 

2.H871E-03 

S.4887E-05 

-6.64HSE-04  3.2086E— 04 

2  d.4b33E-05 

-9.89l5p.-06 

-5.96«7e-06 

-6.37<JlE-06  7.2962E-06 

k 

3  -*.098*t-04 

5.0262F-04 

?.flS97E-04 

1.S803F-05  -4.4021E— OS 

4  -4.H05F.-05 

5.6163F-05 

1 .9120E-06 

1.2830E-07  -2.2846E-06 

j  2,4767F.-03 

-1.4394F-03 

-3.2255E-04 

9.9i«lF-04  -3.8889E-04 

2  -1.0124E-05 

3.4296F-06 

1  .9644E-06 

5.3970F-06  -3.9B24E-06 

3  f 

3.505JF-U4 

-S.60R2F-04 

-8.8412F-05  1.1566E-04 

4  1.84lbE-u5 

-2.H067E-0S 

3.0448E-06 

-1.2226E-07  4.9780E-0? 

(c)  j  = 

=  3 
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Table  XIX 
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Table  XIX  (continued) 


k 


1 

2 

3 

4 

1 

2 

3 

4 

1 

2 

3 

4 


1 

-4.fl751E  00 

-5*9249E-02 

-2.H328F-01 

1.0344E-02 

2.9309K-0? 

-1.6327E-04 

1.164SE-02 

2.1841F-04 

-l.flbbSE-02 

6.2224E-05 

-3.2906F-03 

-8.63H5E-05 


2 

3.8568E-01 

3.5900E-0? 

4.9919E-0? 

-1.6060E-0? 

?.?4fl8E~0? 
-1.0971E-04 
-H.63H7E-0* 
1 #01 17E-05 
-7.2509F-03 

4.2486E-0S 

2. 7208E-03 

-2.  12  79E“0f» 


t 

3 

8.9416E-0? 

-U0355E-02 

-1.1629E-03 

6»3760E-03 

3.30S4E-03 

U5669E-05 

-1.2205E-03 

-2.8250E-05 

8.4902E-04 

-4.6543E-06 

1 .A934E-03 

I.1893E-08 
(b) j  =  2 


4 

-3.H33Mr.-fM 

•1.9107E-4? 

-h.20S7r  -04 

4.4934F-0 J 
-7. 4 3‘JSE-fW 
-H.7ni7fc-0f» 

-4.^4  7SE-04 

4.71  0b£-0h 

l.fl'Vt  -0  < 

2.4943E-llft 

4.20S  U.-04 

-1  ,*>?7M.-48 


5 

-s.*  JS1E-04 

S.21SHE-03 

1  .‘*(10  l*--)4 

-2.3237E-03 

-b. 7‘tPlK-Ob 

-  1  .  ('Mtl1-**  -04 

1 ,*i4M -ob 

o.  j'iO.1  -')4 

2.24il6-oh 
1  .  I  UOM'-llS 

-7.‘SS04E-O  t 
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Table  XIX (concluded) 


1 

-J.4l9ttf-o;» 

-1.6US2E-02 

2  .  fh<ri  \  f — 0  2 

2.64  74F.-0J 

^.4bo2i-:-()3 

"/.flo 5  /t-05 

7  »  4b  1  If— 04 

H.2677E-0S 

-1 .647  7H-0 3 

2.H73UF-o«i 

4,.V>2  If >04 

-3.02806-06 


2 

3.098  7F-03 

1.1232F-02 

4«b6j4F— 03 

-4.4M27E-03 

3.0166E-03 

-S.0343E-06 

"4.41UVI--04 

-fl.b20bf.-06 

S.6‘-M  7F-04 

1  ,fl31 OF-OS 

1.0043F-0J 

6.1M/OF-06 


I 

3 

2.1003E-04 

-3.4617E-03 

3S84F-03 

1.970SE-03 
-?.0«2HE-04 
7. 37236-06 
-1 .693HE-04 

-1 • 092PE-05 

1 .41766-03 

—  3.041 4F  —  06 

?. «3bMt-04 

4.86646-06 


4 

-1 .32276-04 

-2.9578F-03 

4 . 1 Obflf-06 

1.1 712E-03 

-4.6SS5F-04 

-3 .69?3E-06 

2. HO^IF-OS 

2.67S2F-06 

3. K/S4F-04 

S.9468E-07 

-1.22FHF-04 

-1 .1045F-06 


5 

3.8909E-06 

1.S737E-03 

6.97876-06 

— 6.8976E— 04 

-M.Q377E-0F 

-2.1734E-06 

7.0774E-06 

S.4410E-07 

H.4904E-06 

6.5401F-07 

-S.5630t-06 

-2.6257F-07 


Table  XX 

Deflection  Lead  Coefficients 


I  * 

12  3  4 


l 

-3.9351E-02  -1.239HF-0? 

H#9290t-03 

4.68S0F-03 

2 

•S^OOTOE-rtl 

1 .8807F  00 

-2.4089K-02  «U  1390001 

3 

-1.9W?Hh-0? 

7.*  7K7F-0S 

6.721 U “03  -1.S203F-04 

4 

6**»6rtnt>0'3  -3.V^uE-n? 

“1.34S7E-02 

7.4622003 

1 

2.2soot.-o4 

-4.7970F-04 

1  *60940  04 

2 

k 

1.599m£-01  -3.2743F-03 

7.2 O55E-0?  -1.83830  02 

3 

-1 .n?9?t  -»)  / 

‘..l,V‘»f-(>b 

-2.0231F-04 

1.3498  005 

4 

?.684<)F.-(>?  -Y.331?t-'»3 

b.l992F-0j  -3. 9343005 

1 

3.SXS  ft-  — 'tS 

2.2* 78004  — H.4HS0OGS 

2 

-S.  396?f.-02 

1.622HE-03 

-4.U94 3fc-02 

1.0932  0  02 

3 

7.‘V^Mt  -(if  -v,  3  J  JmF-OS 

/.«*U57f«5  -2.0*9hO06 

4 

-1.69311 

S.2n9ot-0.3 

3.«»  13*003  -2.7*(iw003 

(a)  i  =  1 
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Table  XX  (continued) 


i  2 

1  -4,684lE-02  -l#0797F-fl2 

2  2.4711E-02  I.1543E  00 

3  -?.3969E-02  ?.RftQ9F-03 

4  1.6340E-02  -2.985SF-02 

1  ?.7947t-fl4  -9.6383F-05 

2  l,015«E-02  8.4UIE-03 

3  1.6SB6E-06  1.919?F-0ft 

4  l*0662E-03  -2.4102E-03 

1  -1.19361-04  4.6S71F-05 

2  '-3.3991E-03  -5.1704F-03 

3  •P.4R04E-07  8.11M7F-07 

4  -4.702QE-04  1.62S0E-03 

(b) 


3  4 

i.074?E-0?  4.4299E-03  ' 

-3.014RE-02  -8.1854E-02 

i  =  1 

B.l?*0f-03  -l.t39fet-03 
«5*53?3£-03  5.8674E-03 

4 

-5«8<)9iF-«»4  i .68861-04 
•9.3709E-04  -*>.74\«E-03 

i  =  2 

-2.*»V2  7f-04 

J,4989E-03  -}.660lF-03 
/?.8o64t-04  -I.0206E-04  ' 
•1.090IE-03  6.89??£-03 

i  =  3 

l  •  022  1 E-04  -?.6n8F>f-0b 
-4.08491-04  9,73721-04 

4 

=  2 
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Table  XX  (concluded) 


I 

12  3  4 

1  -1.3ttOSfc«*»?  -8.3727E-04  3.P319E-03  6.?4Mf-04 

2  7.919SE-02  1.7003E-02  -P.183AE-02  -4*A02ftF*03 

i  =  I 

3  -7.l96«h-i>3  ?.S07hE-03  ?.4S*3F.-03  -9.0335F-04 

4  3.6638F-03  -3.9977E-03  -].2?3RF-03  1.7?01E*03 

j.  9.o0^0t-ob  “3»4l  1  3f“0S  ”1  •  7808E“04  3*A4bfiE-0*> 

2  4,209«5t-03  6.3B08E-04  -?.41B9E-o3  -2.R182F-04 

k 

3  #j.l»3*fc-«»7  -3.461 3E-0*>  «n.*9SlE-nb  3.?H38E-0b 

4  1.7608F.-04  -3*^7USF-04  2.7399E-04  |'A422F.«ftS 

1  -0*7  1.657lf.»oS  M.R9R*F-ob  -3.06|3f-fl*> 

2  -l.fl31*>L-03  3.4524F-0S  1.9I8SE-03  1.178AF-04 

i  *  3 

3  -3.S444*--ri  f  1.7!]  9F-n*i  U43HF-0S  I768F-0* 

4  O.4209 3.l24bF-04  -?.5?97E-04  •S.122RF-05, 

<c)  j  =  3 
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